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Power Reactor Technology is a quarterly review of developments in reactor technology — 
including new concepts and applications as well as research results, analysis, and experience 
with existing or new components and systems. Prepared by Argonne National Laboratory at the 
request of the Division of Technical Information of the U. S. Atomic Energy Commission, 
Power Reactor Technology is intended for reactor designers and other technical specialists 
concerned with the interplay among reactor development, design, construction, operation, and 
economics. Thus this quarterly journal reports and interprets progress in the power-reactor 
field in terms of its significance to the reactor designer. Its scope includes research and de- 
velopment in the various technical specialties in the field as well as design and construction 
practices, reactor concepts, applications, economics, and operating experience. 

Some articles summarize and critically evaluate reported developments; others review a 
topic broadly so that designers can interpret new developments, experience, and trends. All, 
however, call the reader’s attention to reports and publications that merit study. Because any 
appraisal involves the reviewer’s opinion of the significance of the work reported, readers are 
urged to consult the references for additional information and the judgments of the original 
authors. If a reader has information that causes his evaluation to reinforce, modify, or contra- 
dict the opinions of our reviewers, he is encouraged to write to the Editor. 

For reasonably timely coverage, articles often must be based on fragmentary information. 
Thus, whenever appropriate, continuity will be maintained from one issue to another by use of 
subject headings that reflect parts of the following scope of the journal: 

CONCEPTS and APPLICATIONS: Progress in evaluating the applicability and economics of 
various reactor types and systems (including unconventional approaches), as well as of fuel re- 
sources and cycles, for utility central-station generation ofelectricity, auxiliary power, process 
radiation and heat, desalting, and propulsion—and for terrestrial, undersea, aerospace, and 
other advanced uses. 

ANALYSIS and EXPERIMENTATION: Advancements in the techniques of reactor physics, fluid 


and thermal technology, power conversion, fuel elements, materials, mechanics, control and 
dynamics, reactor and radiation safety, containment, and siting. 


SYSTEMS and COMPONENTS: Experience as reflected in design and construction practice, 
components, systems technology, and operating experience of various specific types of 
reactors —including pressurized- and boiling-water reactors, molten-salt, organic-, gas-, and 
liquid-metal-cooled reactors, as well as generally applicable aspects of research and test 
reactors. 


At Argonne National Laboratory, the Reactor Engineering Division (L. J. Koch, Director) has 
the principal responsibility for the preparation of Power Reactor Technology, with the regular 
assistance of the Reactor Physics Division (R. A. Avery, Director) and the occasional assis- 
tance of other Laboratory divisions (Metallurgy, Chemical Engineering, Chemistry, Reactor 
Operations, and Idaho). Thus, unless noted otherwise, the reviewers and authors of articles are 
Laboratory staff members. The editors will welcome interpretive review articles contributed 
by readers outside the Laboratory. 


JAMES J. DUTTON, Editor; DAVID H. LENNOX, Technical Editor; 
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Ocean Technology and Nuclear Power 


By Edward J. Croke 


Oceanography originated as a basic science 
devoted to the study of the dynamics of ocean 
behavior on and beneath the surface, the map- 
ping of depths and shorelines, and studies of 
the history and composition of the ocean floor 
and undersea ecological system. The fishing 
industry, national defense, and offshore oil 
industry have given rise to an oceanographic 
technology that is based on the application of 
both conventional techniques of naval architec- 
ture and technological developments from the 
aerospace and atomic-energy programs to the 
mastery of the ocean environment. 

The growth of ocean science, technology, 
and industry is leading to a need for power 
sources of various kinds. The nature of these 
needs suggests that the features of nuclear 
power can be adapted in numerous ways to 
enable the exploitation of ocean resources. 


Oceanic Power Needs 


Ocean technology requires power sources 
ranging from data sensors consuming a few 
milliwatts to undersea industrial equipment 
demanding power by tens of megawatts, In 
general, the potential areas of application: for 
nuclear power can be divided into (1) the 
remote and unmanned long-duration systems, 
such as buoys and beacons, that must operate 
for several years without attention and (2) the 
manned and relatively short-duration systems 
that include submersible vehicles and indus- 
trial facilities. The long-duration systems in- 
volve low-power sources in the 1-mw to 100- 
watt range, but the short-duration systems 
tend to have kilowatt or megawatt requirements. 


Thus both isotope and reactor power supplies 
offer promise for many of the projected tasks 
that will be encountered in the conquest of 
inner space. 

Isotepic power sources, such as certain of 
the Systems for Nuclear Auxiliary Power 
(SNAP) systems, can be expected to yield up 
to 100 watts (electrical) or a few thermal 
kilowatts over extended times, whereas isotopic 
steam engines may supply power in quantities 
on the order of 100 kw(t). Nuclear reactors 
could provide power in the range up to tens of 
thermal megawatts. The SNAP-10A (Fig. 1) 
power generator and the naval submarine propul- 
sion plants (Fig. 2) are representative of these 
classes of reactor power sources. 

The potential applications of undersea nuclear 
power are many and varied. Undersea con- 
struction, submarine laboratories, oil drilling 
and production operations, mining and dredging, 
submersible-vehicle propulsion, navigational 
aids, and undersea pumping facilities for pe- 
troleum or water are among the more obvious 
and immediate uses for reliable electric- or 
thermal-power supplies. Thus the future of 
nuclear-powered commercial industry under the 
sea lies in three basic types ofapplications: (1) 
the extractive industries, (2) the harvesting of 
protein resources, and (3) underwater construc- 
tion. 


THE EXTRACTIVE INDUSTRIES 


These include petroleum and mining (ex- 
ploration and production) operations, ocean 
versions of which are becoming increasingly 
valuable. 
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1. NaK outlet from reactor 10. Wires from thermoelectric converter 
2. Thermoelectric NaK pump 11. Batteries 

3. Nak distributor 12. Power outlet 

4. Nak distributor 13. Reactor containment sphere 

5. NaK collector 14. Lead shutdown shield 

6. NaK return 15. Instrumentation and controls sphere 
7. Control drum 16. Lithium hydride 

8. Control drive rod 17. Thermoelectric converter assemblies 


9. Control drive mechanism 


Fig. 1 The SNAP-10A reactor power generator, 
adapted to marine use, is designed to produce 500 
watts with three thermoelectric converter assemblies 
and 2000 watts with 12 converter assemblies; design 
life would be more than five years. 





Fig. 2 Triton, largest U. S. nuclear submarine, is 
propelled by twin reactors of pressurized-water type 
designed for AEC by Knolls Atomic Power Labora- 
tory. 
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Petroleum. Offshore oil production in the 
United States began in 1938 from the Creole 
field off the Louisiana shoreline. The develop- 
ment of the California and Texas offshore 
fields began shortly thereafter. Today, through- 
out the world, underwater oil and gas fields 
are also producing in many inland bodies of 
water, such as Lake Erie, and large-scale 
penetrations of the continental shelves have 
been made in the Persian Gulf, off the coasts 
of Saudi Arabia, in the Gulf of Suez, off Vene- 
zuela, Japan, Mexico, Trinidad, and Peru. As 
extensive exploration of the continental shelf 
areas continues, petroleum production opera- 
tions are expected to expand both geographically 
and in depth. In the essentially unexplored 
coastal depths between 100 and 600 ft, several 
hundred new oil and gas fields are expected to 
be discovered eventually. ! 

Drilling activities at sea have always been 
conducted from fixed offshore towers or from 
floating rigs. These have the problem of station- 
keeping in the face of surface turbulence and 
weather. In the future, undersea drilling from 
submerged platforms or submersibles is en- 
visioned. The freedom of the undersea drill 
rig from the influence of weather, the potential 
for sinking wells at much greater depths (deeper 
than 600 ft) than has been possible with floating 
surfaces, and the resultant extension ofthe area 
available for exploitation all make the prospects 
of undersea drilling most attractive. The only 
“if” is that such operations be proven to be 
economically feasible, which is an almost direct 
function of the technology available. 

A submerged drilling rig with life-support 
facilities has been estimated to require ap- 
proximately 25% more power than acomparable 
surface platform.’ The larger surface rigs 
operating now require about 6000 shp (shaft 
horsepower), so an undersea operation would 
involve a continuous demand of about 7500 shp 
per rig. An additional 500 hp would be required 
approximately every three years for major 
on-site equipment overhauls, and about 10 hp 
would be required for monthly routine main- 
tenance (such as paraffin removal). 


Pressurization of the oil field would be re- 
quired during production. If Seawater were 
used for this, the power requirement would be 
500 to 600 shp; gas pressurization would re- 
quire 2000 to 15,000 shp, depending on the size 
and nature of the field. If the petroleum product 
were transported to shore by pipeline, a pump- 
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ing power ot about 3000 hp would be required, 
depending on local conditions, depth, and pump- 
ing distance. 

The first application of nuclear power to 
offshore petroleum production went into opera- 
tion on June 21, 1965, on anoiland gas platform 
in the Gulf of Mexico. A 60-watt SNAP-7F 
isotopic electricity generator built by the Martin 
Co. (Fig. 3) replaced a diesel-powered gen- 
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Fig. 3 The SNAP-7F isotopic power generator, in 
one of four tubular piles that support unmanned oil 
and gas platform in the Gulf of Mexico, produces 60 
watts to power navigation aids. Design life is five 
years without refueling. 
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erator on an offshore installation owned jointly 
by the Phillips Petroleum and Kerr-McGee 
Companies. As a result of the joint, no-cost, 
research and development contract between the 
AEC and Phillips, the SNAP unit is powering 
two flashing-light beacons and a foghorn during 
a two-year test. This modest beginning might 
be regarded as a first step in the development 
of a nuclear-powered deep-sea petroleum 
industry. 


Mining. Since the British oceanographic ves- 
sel “Challenger” dredged manganese dioxide 
nodules from the sea floor in 1870, abundant 
evidence of the mineral wealth of the oceans has 
been revealed. At present approximately 
$200,000 worth of diamonds are recovered 
monthly along the South African coast; tin is 
dredged from the shallow seas in the Indonesian 
Archipelago; Japan mines iron from its coastal 
waters; and heavy minerals are taken from 
beaches and near-beach areas of the United 
States, Australia, and India. Sulfur is found in 
the Gulf of Mexico, and coal is mined from 
tunnels under the seas off Canada and England. 
Bromine, magnesium, iodine, and common salt 
are recovered commercially from seawater. 
Furthermore, manganese dioxide nodules, phos- 
phorites, gold, platinum, and tin are found on 
the continental shelves and on the deep-sea 
floor. 

The major obstacle to the large-scale ex- 
ploitation of the minerals of the oceans has 
been economic in nature. For example, in 
most cases, the cost of extracting minerals 
from the sea is prohibitively high. All marine 
mining operations have been carried on in 
relatively shallow waters, less than 400 ft deep, 
to allow the use of conventional dredging equip- 
ment.! However, private industry is supporting 
efforts to develop practical methods for under- 
sea mining. 

In a unique arrangement that reverses the 
usual funding pattern, three major firms have 
provided the Bureau of Mines of the Department 
of the Interior with funds and manpower to 
support research and development in marine 
mining systems. The Bureau of Mines has 
supplied project leadership and facilities in a 
research program with the Merritt-Chapman 
& Scott Corp., the Lockheed Missiles and Space 
Company, and the International Mining and 
Chemical Corporation.’ In September 1965, the 
Bureau signed a contract for a cooperative 
exploration of the Pacific continental shelf with 
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the latter two firms.‘ E. P. Wheaton of Lockheed 
estimates that undersea mining stands today 
where the offshore petroleum industry stood 
15 years ago. 


Mineral Exploration. The most immediate 
interest in ocean mining is in alluvial deposits 
close to or within 200 ft of the ocean floor and 
orthogenic deposits (phosphates and manganese) 
on the ocean floor. Nuclear techniques of 
mineral detection through the liquid overburden 
appear both promising and necessary because 
conventional in situ sampling systems cannot 
differentiate between nodular mineral deposits 
and the surrounding sand or mud. In situ 
neutron-activation-analysis devices andneutron 
or beta-particle backscatter detectors must be 
developed to detect and identify mineral de- 
posits qualitatively and quantitatively. Prelimi- 
nary investigation of such detection systems 
is being performed by the Bureau of Mines. 
A nuclear density probe for in-place measure- 
ment of deep-sea sediments has been described.” 

This probe was developed for in situ mea- 
surement of bulk sediment density on the 
deep-sea floor. Measurement is based on the 
scattering of gamma rays from a source in 
the probe. The intensity of the backscatter is 
detected by a G-M counter assembly and re- 
corded photographically by a camera in the 
instrument. The amount of backscatter indicates 
sediment density. After the probe has pene- 
trated the sea floor, the source and detector 
automatically move vertically within the lower 
barrel of the probe and thereby give a density 
profile. Approximately 60 stations can be oc- 
cupied during a single lowering of the device. 
The self-contained probe is approximately 26 ft 
long and weighs approximately 700 lb in air. 

The nuclear technique used in thisinstrument 
involves total absorption or scattering (partial 
absorption) of gamma rays and the subsequent 
detection of the backscattered rays by a detector. 
The source is 100 mc of '*’Cs and the detector 
consists of a bundle of 30 G-M tubes. An 
11.11-in. tungsten—aluminum shield separates 
the source and detector. This device, developed 
by the Lane-Wells Co., under contract to the 
AEC, has proven to be capable of measuring 
densities ranging from 1 to 2 g/cm* with an 
accuracy of 1%. 


Mineral Production. Because open-ocean 
mining involves navigation and site marking, 
undersea beacons are needed to signal the loca- 
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tion of abandoned or temporarily inactive sites, 
to mark the boundaries of mineral deposits, 
and to guide automatic mining equipment. Iso- 
tope-powered sonic buoys that can operate 
continuously for years in the marine environ- 
ment are appropriate for these tasks. Present 
SNAP systems appear overly elaborate and too 
expensive for the purpose; however, they consti- 
tute the best available power sources for the 
durations required. 

Recovery of minerals from the sea floor will 
involve any of several conventional types of 
dredge, suction, or airlift systems. In special 
cases, submersible barges will be needed to 
recover minerals on or beneath the sea floor. 
Water-lift mining techniques are suitable at 
depths less than 1000 ft where power supplies 
in the range of 250 to 1000 kw(e) would be 
appropriate. For deep-water recovery at depths 
from 1000 to 20,000 ft, power sources of 11 to 
30 Mw(e) would be required. ® 

Strip, drift, slope, or shaft mining of the sea 
floor may involve the establishment of semi- 
permanent facilities and perhaps processing 
plants on the ocean floor. An underwater strip- 
mining operation might require 3500 shp for a 
unit capable of handling 900 cu yd of sand or 
similar matérial per hour. Core boresor large- 
scale sample collection would require some 
200 shp for a borer and 5000 shp for a large 
dredger. The processing and cleaning of ores 
would demand power in excess of several 
hundred shaft horsepower; thus undersea mine 
operations can be expected to involve a con- 
tinuous demand for 200 to 5000 shp.” 

For shallow-water mining operations, nuclear 
power sources cannot, at present, compete with 
fossil fuels. For deep-water operations at a 
range of 2100 miles from fossil-fuel ports, 
nuclear power approaches the economic break- 
even point as compared to fossil fuels.° The 
problems of station-keeping, power transmis- 
sion from the surface, surface turbulence, re- 
quired mobility, and power duration may enable 
nuclear power to become competitive with 
fossil fuels on the basis of sheer necessity. 


PROTEIN HARVESTING 


The sea produces approximately 40 million 
tons of protein annually. This yield is being 
harvested by the world fleet oftrawlers, factory 
ships, and whaling vessels that use methods 
analogous to those employed by primitive 
hunters and food gatherers on land. For the 
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fishing industry to improve its efficiency of 
harvest and the yield of sea protein (as it must 
to accommodate an expanding world population), 
scientific “farming” of the sea, aided by modern 
technology, must be implemented. The most 
immediate and promising application of nuclear 
power in the fishing industry is in the supply of 
propulsion and facility power for mother ships 
or factory ships. An estimate based on up- 
grading the largest class of modern mother 
ship indicates that 10,000 shp would be re- 
quired for propulsion and another 10,000 shp 
for electric fishing equipment, cannery power, 
fish pumps, fish-meal processors, and irradia- 
tion equipment.’ 

Small submersibles having crews of two or 
more might be used with the mother ship for 
herding fish or for exploring for biological 
resources. Such craft would have mission dura- 
tions ranging from one day to three months and 
could involve power requirements from 20 to 
several thousand shaft horsepower. Seagoing 
fish-processing plants with refrigeration facili- 
ties, packing capability, or chemical processing 
plants could require 2000 to 30,000 shp, de- 
pending on vessel size, range, and mission.’ 


CONSTRUCTION 


The power requirements for undersea con- 
struction of bridges, piers, or hydrophonic 
platforms are comparable to those associated 
with the extractive industries but require 
shorter durations. 


A VARIETY OF OCEAN POWER NEEDS 


Nuclear instruments such as sediment density 
probes would use an isotopic gamma-ray source, 
such as "Cs, in amounts equivalent to a few 
hundred millicuries. However, sensing instru- 
ments that are not themselves nuclear, such as 
magnetometers, tide gauges, temperature sen- 
sors, seismometers, tsunami detectors, ane- 
mometers, various transducers, and their as- 
sociated amplifiers or telemetry systems, may 
require from a few milliwatts to several watts 
of power over times ranging from a few days 
or weeks in manned research facilities to five 
years for unmanned remote stations. 

Long-term power sources are required for 
unmanned buoys and floating or submerged 
stations that function in such diverse capacities 
as. acoustic beacons for navigation or site 
marking, lighthouses, foghorns, weather sta- 
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tions, communication and surveillance systems, 
recoverable sonobuoys, offshore oil equipment, 
and cathodic protection.® These devices will 
require power from 10 to several hundred 
electrical watts, depending on their specific 
task, and will be required to function without 
attention for up to five years. 

Small manned undersea dwellings for re- 
search or industrial operations will require 
power for heating, cooking, and air regenera- 
tion. Power supplies yielding anything from a 
few kilowatts up to 100 kw(t) are appropriate 
for these purposes. Fortunately much of the 
power required by these stations can be used at 
relatively high efficiency in thermal form rather 
than through thermoelectric conversion sys- 
tems. Heated scuba gear for prolonged work in 
cold water appears to be a possiblility. 

Steam engines in the 1-kw(t) range can be 
used for industrial or research operations 
involving undersea drilling, pumping, valve 
controls, and work-assist devices. 

Small research submarines or industrial 
work vehicles having 1 to 15 shp will be needed 
for manned or remote operation. Seawater 
shielding of nuclear plants offers a possibility 
of obtaining reasonably good power-to-weight 
ratios with such vehicles. Operating times from 
a few days to several months would be required, 
depending on system use. Table 1 summarizes 
these power requirements and indicates the ap- 
proximate duty cycle of potential undersea 
nuclear energy sources. 


Table 1 TYPICAL POWER AND ENDURANCE 
REQUIREMENTS FOR OCEAN POWER? 








Activity Power demand Endurance 
Petroleum 
Drilling 6000 hp + 25% 1 to 3 months 
Workover 10 to 500 hp 1 day to 2 weeks 
Repressurization 500 to 15,000 hp Continuous 
Tanker transport 20,000 hp 1 week/trip 
Pipelining Pumping, depending Continuous 
on design 
Mining 
Dredging 3500 hp 30 days 
Boring 200 to 2000 hp 30 days 
Protein 
Small research vehicles 20 hp 1 day 
Large research vehicles 2000 to 3000 hp 3 months 
Processing vessels 2000 to 30,000 shp 7 to 45 days 
Sea habitats 
Hotel load 10 kw/man 1 to 6 weeks 
Work load Depends on activities 1 to 6 weeks 





Undersea Work Methods 


The extension of man’s ability to work in the 
sea involves the use of four basic methods of 
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operation: (1) the diver, (2) the remotely con- 
trolled and unmanned submersible, (3) the 
manned submarine, and (4) the undersea sta- 
tion. The power requirements of each of these 
can be estimated with some confidence. 


DIVERS 


The diver in a scuba suit or hard-hat rig is 
the most versatile undersea operator in situa- 
tions where a high degree of manipulatory skill 
is involved and risks are acceptably small. 
Recent experience with the Sealab has empha- 
sized the significant problem of maintaining 
body warmth when a man is exposed to the 
ambient ocean; in the low temperatures of deep 
water, the working time of a scuba diver is 
limited to about 30 min. A suit heater that 
would permit a man to remain in the water 
indefinitely would consume about 300 watts 
(thermal).” 


REMOTE-CONTROL VEHICLES 


Remotely controlled submersibles, equipped 
with television cameras and manipulators and 
capable of operating at depths and carrying 
loads greater than a diver could tolerate, can 
be expected to be used in situations that do not 
require a high degree of manipulatory skill or 
in which a diver or manned vehicle would not 
be risked. Examples might be the transport 
and placement of explosives in mining opera- 
tions or the preliminary exploration of an 
undersea site. The robot might be controlled 
through an umbilical cable from a surface 
vessel. If the mission power demand were suf- 
ficiently great, the vehicle could contain an 
integral power supply; otherwise the umbilical 
connection could transmit both control signals 
and power. If the robot contained an integral 
power supply, the umbilical connection might 
be omitted in favor of preprogrammed control 
or sonic control from the surface. The power 
requirement for remotely controlled submers- 
ibles depends greatly on the mission, but a 
fairly large vehicle might require upto 50 kw(e). 


MANNED SUBMERSIBLES 


Submarines for exploration, research, or 
industrial work will probably constitute the 
primary means of exploiting undersea re- 
sources. These vehicles will range in size 
from small one- or two-man vessels with 
mission durations of a few days to large 
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industrial work platforms equipped with cranes 
and manipulators or cargo-carrying submarines 
weighing many tons. Small submersibles with 
crews of one to six persons will cruise at 
maximum speeds of about 8 knots and will 
consume approximately 50 kw(e) for propulsion 
and 10 kw/(e) for life-support facilities. Larger 
vehicles for transporting petroleum or other 
products to shore-based processing plants would 
require 3000 to 5000 shp.? 


Vehicle Design. Although the technology of 
deep-submergence-vehicle design is still rela- 
tively young, more than 22 such vehicles have 
been or are in the process of being built 
(Table 2); Fig. 4 illustrates a typical vessel 
intended for research and exploration. All, 
except the U. S. Navy NR-1, are powered by 
conventional electrochemical sources and are 
limited to a marked extent by the endurance 
capability of their power supplies. For example, 
the bathyscaphe Trieste can remain submerged 
for a maximum of about two days if it does not 
engage in significant maneuvering, or it can 
travel at approximately 1.5 knots for 5 hr using 
a conventional electric-power supply. The first 
deep-submergence research vehicle designed 
to overcome the power-supply limitation will 
be the nuclear-powered NR-1. 


Fisheries Research. A nuclear submersible 
conceptual design was completed recently by 
the Electric Boat Company under the auspices 
of the Bureau of Fisheries Biological Labora- 
tory in Honolulu.'! This agency is charged with 
exploring and developing the fishery resources 
of the central Pacific Ocean. As presently 
understood these resources consist of several 
species of tuna and other large pelagic fishes. 

It has become apparent that neither the tuna 
resource nor the oceanic environment can be 
understood fully from observations made from 
surface vessels. A fishery and oceanography 
submarine was proposed as a better platform 
for tuna research. The known behavior and 
distribution of tuna indicate that such a vehicle 
should have a 20-knot speed, a 1000-ft opera- 
tional depth, and a submerged endurance of 
six weeks. Provisions for direct viewing of 
fish, continuous environmental sampling, and 
quiet operation are essential. 

The feasibility and conceptual design study of 
the vessel shows it to be 163.5 ft long, 23 ft in 
diameter, nuclear powered, manned by a crew 
of 24, and able to carry 7 scientists. Innova- 
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Table 2 EXISTING OR PROJECTED DEEP-SUBMERGENCE VEHICLES"? 
Propulsion 
Operating, Propulsion 
Vehicle name depth, Seed, kaa endurance, No. of propulsion Energy aes 
Owner/operator ft Cruise Max. hr/at knots units source Length Beam Height 
Aluminaut 
Reynolds International, 15,000 3 3.8 32/3 Two 5 hp horiz. (tail Four 77-cell 55 8 14.25 
Inc. (1965) plane); one 5 hp vert. 32,500 amp-hr 
(amidships) silver— zinc 
(internal) 
Alvin I 
Office of Naval Re- 6,000 2.5 Q 8/2 Two 5-hp hydraulic Three external 22 8 13 
search/Woods Hole motors powered by oil-encased 
Oceanographic pumps driven by two lead-acid 
Institution (1965) oil-encased electric pressure 
motors compartments 
American Submarine Co. Model 300 
American Submarine 300 1 6 3/3 One 3 hp; 12 to 24 Four 6-volt 13 4.2 4.75 
Co. (1961) volts lead-acid 
American Submarine Co. Model 600 
American Submarine 600 1 6 6/3 One 3.5-hp oil- Four 6-volt 13 5.5 5.2 
Co. (1963) immersed motor nickel— cadmium 
Archimede 
French Navy (1961) 36,000 3 One 30 hp, 110 volts Two 110-volt 860- 69 13 26.5 
(horiz.); one 5 hp amp-hr oil- 
(vert.); one 6 hp encased alkaline 
(direction control) batteries (out- 
side pressure 
hull) 
Auguste Piccard 
Swiss Natl. Expos. 2,500 lto2 6 8/6 One 80 hp Two 2000-amp- 93.5 19.7 24 
Corp. (1964) hr 110-volt lead— 
acid batteries 
Benthos V 
Lear Siegler, Inc. 600 2 2/2 Two 1 hp, 24 volts One silver— 11.3 6.1 6 
(1963) cadmium 
Cubmarine PC3A ’ 
Perry Submarine 300 3 4.25 8/2; One 7 hp Six 6-volt 210- 18.5 3.5 6 
Bldrs., Inc./Ocean 3/3 amp-hr lead—- 
Systems, Inc. (1964) acid 
Cubmarine PC3B 
Perry Submarine 600 2 4 2/4; One 7 hp, 36 volts Fifty 1.5-volt 22 3.5 6 
Bldrs., Inc./Ocean 4/2 250-amp-hr 
Systems, Inc. (1964) silver— zinc 
Deep Jeep 
Naval Ordnance Test 2,000 2 2 4-6/2 Two 0.75 hp, 24 volts Eight 6900-watt- 10 8.5 8 
Station, China Lake, (variable); oil- hr lead-acid 
Calif. (1964) immersed 
Deep Star 4000 
Westinghouse Electric 4,000 1 3 6/3; Two 4.5 hp a-c free- Three lead-acid 18 11.5 7 
Corp. (under con- 12/1 flooded 400-amp-hr in 
struction) 62 2-volt cells 
(inverter) 
Diving Saucer 
Westinghouse Electric 1,000 1 1 4/1 One 2-hp motor driving Three 104-amp- 9.5 9.5 5.5 
Corp. (1959) centrifugal water hr lead-acid 
pump; two jets pro- batteries 
vide thrust 
Dolphin AG(SS)555 
U. S. Navy (under con- One diesel-electric Three silver— 152 19 
struction) zinc batteries; 
two diesel— 
electric 
Kuroshio II 
Hokkaido Univ., Japan 650 2 Tethered One 440-volt a-c Electric from sur- 36.7 7.15 10.4 
(1960) motor face through 
cable 
Moray 
U. S. Naval Ordnance 6,000 6 15 3.6/6 One 90-hp battery- Bank of 240 33 5.3 5.3 
Test Station, China operated torpedo silver— zinc 
Lake, Calif. (1964) motor; counter-ro- secondary cells 
tating propellers 
Severyanka 
USSR/All Union Inst. 550 15 16,500-mile Diesel—electric Diesel— snorkel; 240 22 15 
of Marine Fishery range lead—acid 
and Oceanography 
(1958) 
Star I 
General Dynamics 200 0.75 3/0.75 Two side-mounted; Two external 18- 10.1 6 5.8 
Corp. (1963) 0.25 hp, 18 volts volt d-c lead— 
acid 
Star II (Asherah) 
General Dynamics 600 1 4 10/1; Two side- mounted; External 24 volt 17 7.7 7.6 
Corp. (1964) 2/4 2 hp, 24 volts lead-acid 
Star III 
General Dynamics 2,000 1 4 10/1; 7.5-hp single-screw 115-volt 30-kw-hr 24.5 6 8 
Corp. (under con- 2/4 stern drive; one oil-immersed 
struction 2-hp bow thruster; lead-acid 
one 2-hp vert. 
hovering motor 
Trieste II 
U. S. Navy (1964) 20,000* 2 5/2 Two 10-hp aft; one External 145 kw- 67 15 18 
2-hp bow thruster hr lead-acid 
Yomiuri 
Yomiuri Shimbun 1,000 4 6/4 One electric (a-c); Diesel-electric; 48 8.2 9.2 
Newspaper, Tokyo, two diesel auxiliary diesel 
Japan (1964) for battery charg- 
ing; alternator for 
a-c power 
NR-1 
U. S. Navy Continental Indefinite Two or more Pressurized-water 
shelves nuclear reactor 





*With Terni pressure sphere; 36,000 ft with alternate Krupp sphere. 











Fig. 4 Research submarine designed to take 7000-lb 
payload to 6000-ft depths. This 50-ton 40-ft-long 
electrochemically powered vehicle is being developed 
and built by Lockheed Missiles and Space Company. 


tions include an observation sphere, 2 research 
laboratories, several deep-sampling devices, 
and a 4-pod motor arrangement. Table 3 sum- 
marizes the characteristics of this vehicle. 
Argonne National Laboratory (ANL) also 
studied a nuclear submarine concept for ocean- 
ographic research.” The ANL design featured 


Table 3 PROPOSED FISHERIES-RESEARCH 
NUCLEAR SUBMARINE"! 





Hull material HY-80 steel 
Hull diameter 23.0 ft 
Length, overall 163.5 ft 
Power 3240 shp total 
Per motor pod 810 shp 
Reactor 17.0 Mwit) 
Displacement, normal surface 1265 tons 
Net submerged 1407 tons 


Stability and trim lead 
Margin lead 


37 tons (2.9%) 
145 tons (11.5%) 


Main ballast 142 tons 

Reserve buoyancy 11.2% 
Maximum speed, submerged 20 knots 

Surfaced 11 knots 
Turning radius at slow speeds 150 ft 
Depth, operating 1000 ft 

Collapse 1500 ft 
Dive duration 9 weeks 
Forward laboratory space (cylindrical) 

Volume 2900 cu ft 

Floor area 290 sq ft 
Aft laboratory space (cylindrical) 

Volume 1450 cu ft 

Floor area 65 sq ft 





a 7.75-Mw pressurized-water reactor (Fig. 5) 
that could supply 2000 shp for the propulsion 
cf a 100-ft-long 15-ft-diameter vehicle. 


UNDERSEA STATIONS 


Small undersea facilities, such as ocean- 
floor laboratories, will require thermal and 
electric power for lighting, heating, cooking, 
and life support. The total power demand is 
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1200 psia, 2000 shp 
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3.2 psia, 
145°F, 
14% moisture 
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Fig. 5 Fisheries-research submarine concept stud- 
ied by Argonne would be propelled by 7.75-Mw pres- 
surized-water-reactor system; shield diameter would 
be 10 ft and propulsion system would be 27 ft long 
overall.?? 


expected to be approximately 25 kw(e) for a 
crew of 10 divers. 

Large undersea processing plants for the 
separation of oil and gas or the benefaction 
of ores might require several thousand kilo- 
watts of energy for hotel functions alone.” 
Integral or detached stationary power supplies 
(Fig. 6) will accompany such facilities. 





Fig. 6 Typical undersea utility plant envisioned by 
Aerojet-General Nucleonics might supply electric 
power, potable water, interior heating, hydraulic 
power, and oxygen for undersea plants that process 
oil, gas, oY ores. 


Oceanic Power Sources Available 


Power sources for undersea use are limited 
to methods that neither produce exhaust fumes 
nor require abundant air suppiies but which can 
supply relatively large quantities of power. 
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Thus batteries, fuel cells, and nuclear sources 
are the only reasonable possibilities at this 
time. Both endurance and the power required 
must be considered when the various chemical 
and nuclear sources are compared in the light 
of anticipated oceanic missions. 


Comparison. Figure 7 provides a compari- 
son among various 50-kw(e) power sources in 
terms of system weight and operating time.’ 


1500 — 























Lead-acid 
battery Silver-acid 
2 battery 
= 
=|5 1000 3 
=|8 50-kw/(e) boiling - water reactor 
z 2 50-kw(e) thermoelectric nuclear reactor 
am = 50-kw(e) isotopic thermoelectric 
B15 500K 50-kw/e) fuel-cell plant (without fuel 
= [ containers) 
= 
Y Peas | 
jd 
I | | | | | ! J 
0 20 40 60 80 
Time at Full Power, hr 
Fig. 7 Comparison of system weights and operating 
time of various 50-kw(e) power sources.? 


Thus, at the 50-kw(e) power level, it is evident 
that fuel cells and isotopic electric systems 
offer a clear weight advantage for missions 
longer than one day. 

For higher powers, such as those associated 
with large submersibles, or for undersea tasks, 
such as mining operations, where there is a 
continuous high power demand over an extended 
time and weight is not of major significance, 
nuclear sources are more than competitive. 
Although the power threshold at which the nu- 
clear reactor becomes competitive with other 
sources on a power-to-weight ratio basis has 
not been absolutely established, a report!® of 
the Interagency Committee for Oceanography 
of the Federal Council for Science and Tech- 
nology concludes that, for large undersea ve- 
hicles requiring power plants of 1000 shp or 
more and an endurance greater than 400 hr, 
nuclear power is “indicated.” Even for smaller 
vessels with more modest requirements, 30 shp 
and 30 hr, the report concludes that batteries 
usually will not be practical. The report finds 
that the alternatives are various types of 
power sources, including nuclear sources, all 
of which would require considerable develop- 
ment effort. 
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Nuclear Power in the Sea. Among the means 
considered for supplying power over the broad 
range of undersea technology are a number of 
thermomechanical, thermoelectric, and turbo- 
electric energy systems that use isotope or 
reactor power. Their relative merits depend 
strongly on the weight, duration, shielding, 
and power required for a given mission or 
task. A number of nuclear power supplies 
have already been used in oceanographic mis- 
sions, and other nuclear oceanic systems are 
being developed. 


ISOTOPES AT SEA 

Relatively few isotopes satisfy the criteria 
for use aS marine power sources. These cri- 
teria include a half-life exceeding 100 days, 
a thermal power density greater than 0.1 
watt/g, ready availability, and relatively low 
cost. The AEC is focusing its attention on four 
isotopes that meet these requirements: ‘Sr, 
Nios, M4ce, and Co. Table 4 shows the 
pertinent characteristics of these isotopes. 
Plutonium-238 and ‘’Pm are also possibilities 
for limited application. It appears that the na- 
tional output of these isotopes will be such 
that several thousands of kilowatts of isotope 
power will be made available annually. AEC 
isotope-application programs span a power 
range from 20 puw(e) to 10 kw(e) and include 
thermal, thermoelectric, and thermodynamic 
energy-conversion systems. 


Isotopic Thermal Power. A joint AEC—Air 
Force program to investigate the application 
of isotope power to life-support systems, in- 
cluding water recovery, carbon dioxide re- 
moval, oxygen recovery, heating, cooling, and 
particle removal, is being conducted by Gar- 
rett Corporation. The system is designed to 
use a radioisotopic heat source. Isotopes be- 
ing evaluated for this program are 7**pu at 
Mound Laboratory and '’Pm at the Hanford 
plant. Management of the development effort 
is performed by the Thermal Applications 
Branch of the AEC Division of Isotopes Develop- 
ment (DID). Although this program is primarily 
oriented toward space-flight applications, the 
system also is regarded by the DID as a 
possible undersea-laboratory life-support sys- 
tem.’ 

Evaluation of a proposed isotope-powered 
thermoacoustic system for an undersea sonic 
beacon is in progress at the DID, which is 
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Table 4 ISOTOPIC POWER SOURCE DATA*® 








%Sr-y BICs—Ba ‘M44Ce—Pr 47pm 80C5 

Half-life, years 27.7 30 0.78 2.6 5.2 
Radiation B,y By B By 
Power, curies/watt (thermal) 148 207 126 2790 65 
Typical fuel form SrTiO, Glass CeO, Pm,0,; Metal 

Density 3.7 3.1 6.3 6.6 8.7 

Power, watts (thermal)/em3 0.8 0.24 20 1.8 27 
Possible availability, kw(t)/year 

1965 £2 5 8 0.4 100 

1970 90 4 800 10 150 

1980 650 620 4200 56 
Cost, $/watt (thermal) 

1965 AEC price 112 104 126 4900 33 

1966 AEC price” 30 26 19 556 33 

Future estimate 19 21 0.90 93 10 
Cost, $/kw-hr 

Thermal 0.40 0.40 0.20 15.00 0.40 

Electrical 10.0 10.0 30.0 10.00 





considering a study of candidate systems for 
directly converting heat into sound. If such a 
system proves feasible, it would be less costly 
and complex than comparable heat —electricity— 
acoustic energy systems. A system proposed 
by the Martin Co. envisions the use of heat 
from a Sr source to boil water. The steam 
would power a hydroacoustic transducer that 
produces sound from the flow of oil through 
an orifice in a diaphragm. The entire system 
could fit into a lightweight 1-cu ft box and has 
a projected lite’ of 10 years. 

An isotope-heated scuba suit, being devel- 
oped at Mound Laboratory, combines a 300-watt 
238py heat source with a fluid-flow heat ex- 
changer; the fluid would transmit the heat 
through veinlike tubes embedded in the scuba 
suit. This heating unit will weigh approximately 
15 lb and is designed to be backpacked. It will 
extend the useful immersion time of a diver 
from approximately 30 min/day up to the limit 
of his physical ability to perform underwater 
work.’ This development program also is being 
sponsored by the Thermal Applications Branch 
of the DID. 

A possible application of Sr as an undersea 
heat source is being investigated at Oak Ridge 
National Laboratory (ORNL). Although this 
heater, essentially an encapsulated isotope brick 
pile, has been proposed as a heat source for 
the second phase of the Sealab experiment, 
no decision to develop or use it has been 
announced. 

In view of the prospective requirements for 
larger sources of thermal energy for undersea- 
facility heating and power sources, it is rea- 
sonable to expect that the present relatively 


modest AEC effort along these lines will be 
expanded as the oceanographic (and aerospace) 
programs progress. 


Isotopic Thermoelectric Systems. The AEC 
has concentrated considerable effort on the 
development of isotopic thermoelectric sys- 
tems. A family of such power supplies has 
been built and is being used in various marine 
capacities; more advanced systems are also 
being developed. 

Several prototype "Sr-fueled lead telluride 
thermocouple systems have been placed in op- 
erational test by the AEC. These generators, 
spanning the power range of 7.5 to 60 watts, 
were developed under the SNAP-7 program by 
the Martin Co. 

The SNAP-7A device (Fig. 8), a 10-watt gen- 
erator that supplies power for a Coast Guard 
navigational buoy at Curtis Bay, Md., contains 
40,000 curies of "Sr in the form of strontium 
titanate pellets. This fuel is contained within 
two layers of Hastelloy C to form a total 
thickness of 0.5 in. Sixty pairs of lead tel- 
luride thermocouples, having a hot-junction 
temperature of 800°F and a cold-junction tem- 
perature of 115°F, supply 10 watts at 5 volts 
direct current. The entire unit, including bio- 
logical shielding, weighs 1840 lb, of which the 
fuel, generator, and containment capsules weigh 
only 121 lb. 

The SNAP-7B (Fig. 9) is a 60-watt lead 
telluride thermocouple generator fueled with 
225,000 curies of strontium titanate. The fuel 
is encased in a double layer of Hastelloy C and 
shielded by 3 in. of depleted uranium. The unit, 
including shielding, weighs 4600 lb. The SNAP- 


Winter 1965-1966 


Electric lantern 














26.5 ft 
tte 
— ~~ 
oO 
—yw 
—~“ 

4 

re 

re 











a 
Batteries 
Fe aE 3 
— 


10-watt strontium -90 
thermoelectric generator 








\ 





t 
i‘ if 


N 


Fig. 8 The SNAP-7A isotopic thermoelectric power 
system, in base of buoy, generates 10 watts to re- 
charge batteries in Coast Guard navigational buoy at 
Curtis Bay, Mad. 


7B is in service as the power supply for the 
Baltimore lighthouse in Chesapeake Bay. 

The SNAP-7C system, a 10-watt unit that is 
essentially identical to the SNAP-7A, supplies 
power for a U. S. Navy automatic weather 
station in the Antarctic. 

The SNAP-7D, a 60-watt generator similar to 
SNAP-7B, powers a floating Navy automatic 
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Fig. 9 The SNAP-7B isotopic therimoelectric power 
system has 120 pairs of lead telluride therimocou- 
bles to convert heat from 14 fuel capsules to 60 watts 
that power Coast Guard’s Baltimore Light. 


weather station in the Gulf of Mexico (Fig. 10). 
This buoy system, called NOMAD (Navy Ocean- 
ographic and Meteorological Automatic Device), 
is a forerunner of a network of such installa- 
tions planned by the Navy. 

The SNAP-7E (Fig. 11), which is essentially 
similar to the 7A and 7C units but with an out- 
put of 7.5 watts, provides power for a naviga- 
tional sonobuoy on the 15,000-ft-deep floor of 
the Atlantic Ocean near Bermuda. Unlike the 
7A and 7C units, SNAP-7E has a forged-steel 
pressure vessel to withstand the 8000-psi am- 
bient pressure at its operational site. 

The SNAP-7F generator (Fig. 3) is a 60-watt 
device that has powered navigational aids on an 
offshore oil and gas platform in the Gulf of 
Mexico. This operational test unit is similar 
to the other 60-watt generators in the SNAP-7 
series. 

A second generation of marine SNAP sys- 
tems to produce powers of 10 to 60 watts is 
being developed. These are the SNAP-21 units 
that will operate in pressures up to 10,000 psi 
for at least five years without maintenance. 














Fig. 10 The SNAP-7D isotopic thermoelectric gen- 
erator provides 60 watts to power weather instru- 
ments and telemetry system on a Navy unmanned 
weather boat in the Gulf of Mexico. 


The SNAP-21 units will be developed for 10-, 
20-, and 60-watt powers and will be fueled 
with "Sr. The Minnesota Mining & Manufactur- 
ing Co. has undertaken the development effort, 
which is scheduled to produce a 10-watt proto- 
type unit in 1966. The use of improved ther- 
moelectric materials, higher hot-junction tem- 
perature, and improved thermal insulation and 
design will result in second-generation sys- 
tems whose cost and weight will be an order 
of magnitude less than the SNAP-7 series. 
For example, the 7.5-watt SNAP-7 unit weighs 
6000 lb, but the 10-watt SNAP-21 will weigh 
only about 500 lb (Ref. 15). 

The AEC also is planning an Advanced Large 
Milliwatt Generator Program. The first phase 
o: the program would yield conceptual designs 
for 150-, 600-, and 900-mw generators that 
could power sonic pingers for wellhead locators 
and other such low-power devices. 

The SNAP-15 program involves development 
of a 1- to 3-mw thermoelectric generator 
weighing less than 1 lb. Such a power con- 
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Fig. 11 The SNAP-7E isotopic thermoelectric gen- 
erator, which is on sea floor at 15,000 ft depth in the 
Atlantic Ocean, produces 7.5 watts to power a Navy 
navigational sonobuoy that is about 3000 ft above the 
generator. 


stitutes the lower useful limit for undersea 
energy sources, but SNAP-15 could power a 
variety of instruments and data sensors. The 
General Atomic Division of General Dynamics 
Corp. is the contractor for the development of 
this ***pyu-fueled generator. 


Isotopic Engines. Thermodynamic engines 
that can supply power in the range 10 to 100 
kw(t) will find use in small submersible propul- 
sion systems, turbodrills, and work-assist de- 
vices. 

A conceptual design study of an isotopic heat 
source integrated with a propulsion engine for 
undersea craft is in progress at Aerojet- 
General Nucleonics under a cooperative con- 
tract with the DID. This RISE (Radioisotope 
Submersible Engine) program will provide the 
first demonstration of a thermodynamic isotopic 
ocean power system having a cycle efficiency 
two to four times that of existing thermoelec- 
tric conversion devices.® 

The RISE concept, for which ®°Co has been 
tentatively selected as the heat source, envi- 
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sions a _ turbine-Rankine-cycle system using 
Dowtherm A as the working fluid in a direct 
(single-loop) power-conversion system that will 
consist of a combined rotating unit that includes 
a turbine—alternator, condensate pump, regen- 
erator, and condenser. The power-conversion 
system will be sealed to prevent leakage of 
working fluid or contamination of the system. 
Current effort in the RISE program involves 
the conceptual design, detailed program plan- 
ning, and initial research and development for 
a 3-kw(e) or 4-shp demonstration engine. The 
System design, like the heat source, will be 
scalable to about 20 kw(e) or 26.8 shp. This 
proof-of-principle engine probably will use a 
16-kw(t) isotope source with design potential 
for being upgraded’ to 100 kw(t). 

In an associated effort the Savannah River 
Laboratory is performing a source-selection 
trade-off study, including research and devel- 
opment, on the production and use of °Co for 
this purpose. A study of Sr as an alternative 
heat source is in progress at ORNL. 

An isotope-heated steam piston engine also 
has been studied at Aerojet-General Nucleon- 
ics.1° 17 This system features a 12-shp two- 
cylinder steam engine weighing approximately 
3800 lb, including 2900 lb of shielding. Figure 
12 and Table 5 illustrate the system and 
delineate its characteristics. Isotope steam 
engines Similar to this propulsion power pack- 
age are also being studied for underwater 
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Fig. 12 Radioisotope-powered steam piston engine 
could propel small submarines or power turbodrills 
or other work-assist devices. 
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Table 5 SUMMARY OF RADIOISOTOPE-ENGINE 
CHARACTERISTICS"? 





Cycle characteristics 
Shaft power output 
Maximum pressure 
Maximum temperature 
Minimum temperature 
Ideal efficiency 

Boiler characteristics 
Fuel 
Cladding material 
Fuel-element configuration 
Pin size 


Total fuel volume 
Maximum thermal power 
Minimum thermal power 
Shield material 
Engine characteristics 
Type 
Speed 
Auxiliaries 


Condensers 


Feed pumps 


12 kw 
1200 psia 
567°F 
141°F 
34.8% 


Cerium-144 as CeO, 
Hastelloy C 

19 pins in hexagonal array 
1.5 in. in diameter by 7.9 


in. long 
250 cu in. 
71 kw 
45 kw 


Natural or depleted uranium 


Reciprocating, double acting, 
compound 2-cylinder 
2000 rpm 


Natural convection to sea- 
water 
2 steam-driven reciprocat- 


ing pumps; 2 electrically 
driven gear pumps 





mechanical power supplies and work-assist 
devices. 


REACTOR-POWERED SYSTEMS 


The high power range of reactor systems has 
been effectively exploited for the propulsion of 
ships and large submarines. However, for the 
smaller research submersibles intended for 
deep penetration of the ocean, no suitable 
reactor power or propulsion systems are avail- 
able now. Vehicles ranging in size between the 
large combat submarine and the very small 
(< 300 tons) submersibles would lend them- 
selves well to reactor power. Of interest for 
this application is a class of reactors in the 
SNAP-10A family, which was originally devel- 
oped for space-flight applications and generates 
500 watts (thermal), and the higher powered 
SNAP-8, which produces 100 kw(t) and could 
generate 20 to 50 kw(e). 


SNAP-10A is a uranium—zirconium hydride- 
fueled reactor that contains 4 kg of **U. 
Weighing approximately 1000 lb, SNAP-10A 
measures 5 ft at the base and is 12 ft high. 
Its NaK coolant is circulated by a thermo- 
electric pump, and power is converted by 
silicon—germanium thermoelectric mod- 
ules with 885°F hot junctions and 235°F cold 
junctions. 

An adaptation of SNAP-10A that would yield 
500 to 2000 watts (electrical) has been studied 
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for stationary undersea use (Fig. 1). Total 
weight of the 2000-watt version of this marine 
reactor is 8200 lb (with shielding for retrieval); 
the 500-watt version would weigh 7200 lb. The 
SNAP-10A reactor was developed by Atomics 
International under contract to the AEC. 

Larger versions capable of delivering 300 to 
2000 kw(e) have been developed and might be 
adapted for undersea use. Such reactors could 
be of service in large operations such as 
undersea construction; benthic laboratories; 
offshore oil drilling; and production, mining, 
and dredging. 


The NR-1 manned submersible represents 
the only large-scale effort now being made to 
develop a reactor-powered deep-sea research 
vehicle. Although it is expected that the proven 
nuclear technology of the combat submarines 
will contribute materially to the configuration 
of the NR-1 pressurized-water reactor, little 
power-plant design information is available. 


Supercritical -Water Reactor. A reactor con- 
figuration that might be appropriate for propul- 
sion and auxiliary power inadeep-submergence 
vehicle was studied by Battelle Memorial Insti- 
tute.!® This supercritical-water reactor (SWR) 
plant, whose design is summarized in Table 6, 
has a number of features that make itattractive 
as an undersea power supply. The heat capacity 
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of supercritical water averages twice that of 
ordinary water and six times that of sodium. 
This large heat capacity, coupled with a large 
temperature rise, results in low flow require- 
ments and thus smaller pumps and piping. 

Water is a good neutron moderator and good 
heat-transfer medium; the relatively high fluid 
temperature associated with the system leads 
to high efficiency and a lightweight power- 
conversion system. Moreover, the SWR is a 
natural extension of the pressurized-water- 
reactor concept for which considerable design 
and operating experience exists. The usual 
requirement for a thick pressure vessel for a 
water reactor is less severe in the SWR design 
because ocean pressure can be used topartially 
balance working-fluid pressure at operating 
depths. 


SWR Vehicle. A research vehicle was con- 
sidered for missions to 10,000 ft. Composed 
of three adjacent 10-ft-diameter spheres with 
the supercritical-water-reactor power package 
suspended below them, the entire system would 
be surrounded by a hydrodynamic firing. One 
of the spheres houses the operating crew, 
another is for the research crew, and the 
third is for personnel-support equipment. Ra- 
diation shielding is supplied by approximately 
8 ft of seawater in addition to the reactor 


Table 6 SUPERCRITICAL-WATER REACTOR DESIGN DATA!8 





Core 
Active diameter 
Height 
Structural material 
Core power 
Reflector thickness 
Control 


Fuel elements 
Type 
Meat 
Fuel coupon thickness 
Cladding thickness 
Coolant-channei width 
Element size 
Plates per element 
Number of elements 
Initial fuel loading 
Maximum meat temperature 
Maximum surface temperature 
Heat-transfer data (design conditions) 

Coolant 

Pressure, average 

Inlet temperature 

Outlet temperature 

Core pressure drop 

Flow rate 
Core flow area 
Heat flux (average) 


12 in. equivalent 


18 in. 
Type-347 stainless steel 
1440 kw 


2.5 in. 
Stainless-steel-clad boron carbide 
vanes in reflector 


UO, coupons in stainless-steel plate 
93% dense UOg, fully enriched 

40 mils 

7 mils 

15 mils 

1.125 in. square, including 20-mil can 
15 

45 

64.8 kg of 235y 

1170°F 

1140°F 


4340 psia 

300°F 

1000°F 

<1 psia (without orifices) 
4400 lb/hr 

0.07 sq ft 

29,000 Btu/(hr)(sq ft) 
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thermal shield and pressure vessel. A reactor 
and power-conversion package providing 250 
kw(e) for propulsion and auxiliary electric 
power for extended missions has a mass of 
approximately 15,700 lb at a specific gravity of 
approximately 1.8. 


Surface-Ship Propulsion. Reactor propulsion 
for surface vessels has been a reality for 
some time. The nuclear merchant ship Savannah 
(Fig. 13), the aircraft carrier Enterprise (Fig. 
14), the cruiser Long Beach, the destroyer 
Bainbridge, and the Russian icebreaker Lenin 
have amply demonstrated the operational feasi- 
bility of nuclear surface propulsion. 


Gas-Cooled Reactor. The gas-cooled type of 
reactor has been evaluated by several organiza- 
tions as a promising hydrospace power source 
for stationary installations. An Aerojet-General 
Nucleonics study’® indicates that, for near-term 
and low-power (<1 Mw) requirements, an opti- 
mum configuration would consist of a hetero- 





Fig. 13 The N.S. Savannah, a 22,000-ton vessel, was 
the first cargo-passengey ship to be propelled by nu- 
clear power. 
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Fig. 14 The USS Enterprise has demonstrated bene- 
fits of nuclear power for aircraft carriers. 
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geneous light-water-moderated reactor cooled 
by a breathing-gas mixture (helium, nitrogen, 
and oxygen). In addition to serving as both the 
reactor coolant and as the hull-pressurizing 
medium, such a gas mixture could be replaced 
on site by use of the breathing-gas generation 
system. Power would be converted by a Brayton- 
cycle system with a split-turbine-driven alter- 
nator and recuperator. Plant output would be 
100 kw(e) net plus 100 kw(e) for auxiliaries 
and hotel load. A precooler would provide the 
cycle heat sink; a natural-convection seawater 
heat exchanger would provide the ultimate heat 
sink. The reactor compartment would be flooded 
to provide shielding for the power-conversion 
compartment. 

MH-1A. A possible precursor of the large 
industrial undersea nuclear power plant of the 
future is the Army MH-1A floating nuclear 
power plant (Fig. 15). This 10-Mw/(e)-reactor 
electricity-generating system was designed and 
built by the Martin Co. under contract to the 
Army Corps of Engineers. Installed ina Liberty 
Ship hull, the MH-1A will provide a floating 
power supply that can be towed to any port in 
the world to supply or augment electric power 
for emergency conditions. Table 7 summarizes 
system characteristics.”° 

TURPS. A somewhat smaller reactor ther- 
moelectric power system is being studied by 
the Martin Co. under contract to the Surface 
Nuclear Power Group of the Air Force Weapons 
Laboratory. This TURPS (Terrestrial Unat- 
tended Reactor Power System) was proposed 
as a power source for Air Force remote radar 
installations; however, it also is being evaluated 
as an undersea power plant. 

The TURPS concept has a uranium —zirco- 
nium hydride reactor that is cooled by natural 
circulation of tetraphosphorus trisulfide. Its 
rated thermal power is 1.56 Mw, or 100 kw(e). 
Shielding is accomplished with earth or water 
and two shield plugs above and below the 
reactor; a shield sleeve surrounds the reactor 
during transportation. The cylindrical core 
measures 13 by 13 in., and the entire assembly 
is 4 ft in diameter and 19 ft long. The weight 
of the system”! is 18,700 lb. According to the 
Martin Co., TURPS would be about 200 times 
more powerful than any direct- conversion sys- 
tem now operating and would be able to operate 
for five years at full power without refueling. 

From this enumeration of reactor power 
plants can be seen the considerable variation 
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Fig. 15 The MH-1A floating 
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Table 7 MH-1A FLOATING NUCLEAR POWER PLANT?” 


Length: Overall 441.5 ft, beam 65 ft 

Displacement: 8254 tons 

Gross output: 45 kw(t) 

Net output: 10,000 kw(e) 

Containment: Reactor system is com- 
pletely contained in steel contain- 
ment vessel; additional safety is 
provided by special collision and 
grounding barriers in hull 

Shielding: 600 tons of lead and 
polyethylene and 880 tons of con- 
crete surrounding the steel contain- 
ment vessel 

Plant design life: 20 years 

Core type: stainless-steel tubes con- 
taining UQ, fuel pellets 

Reactor type: pressurized water 

Fuel elements: 3328 stainless-steel 
tubes 39 in. long and 0.507 in. in 
diameter; fuel pellets are low- 
enrichment UO, 

Core size: 30 in. high, 45 in. in equiv- 
alent diameter; 36 in. in active 
height 
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nuclear power plant provides 10 Mw(e) to any area to which a barge 
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Control rods: 12 cruciform cross- 
section rods of boron—stainless- 
steel alloy 
Core life: 1 year at full power before 
reshuffle and partial replacement 
Uranium loading: 1460 kg with 4.65% 
enrichment and 1460 kg with 4.07% 
enrichment 
Primary-loop operating pressure: 
1500 psia 
Average temperature: 
Secondary system 
Steam pressure: 330 psia (throttle) 
Steam generator: Verticle U tube 
Turbine-generator: 3600 rpm 
direct-connected; 13,529 kva at 
60 cycles 

Electrical: 10,000 kw(e) net at 60 
cycles and 13,800 volts; 8083 kw(e) 
net at 50 cycles and 11,500 volts 


190°F 
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among the reactor configurations that are being 
considered for use in the ocean. Except for 
naval reactors, all the systems available now 
are essentially land-based systems modified 
for a maritime setting. Until additional ex- 
perience with nuclear power in an oceanographic 
framework is accumulated, it will be difficult 
to identify the specific characteristics and 
peculiarities of an optimum undersea nuclear 
power plant. 

The AEC has initiated a program for the 
development of low-power nuclear reactor sys- 
tems that initially will cover the 100-kw/(e) to 
2-Mw(e) range and be applicable to a variety of 
uses, such as oceanographic power plants.”? 


Conclusions 


The development of nuclear power sources 
for use in the ocean is progressing steadily 
along a broad, if somewhat discontinuous, 
front. In the extremely low-power range, nu- 
clear or nuclear-powered instrumentation tech- 
nology has received less attention than it 
deserves. Isotope-thermoelectric intermediate- 
power source development is well into the 
proof-of-principle stage and is being advanced 
in a vigorous and systematic way, probably be- 
cause of the impetus of the aerospace program. 

Although isotope thermal and thermodynamic 
power sources have been investigated, experi- 
ments have been performed, and analyses have 
been made, most of these devices are still in 
the laboratory or conceptual-design stage and 
much work remains to be done. The lack of a 
specifically defined immediate requirement for 
such systems has probably delayed their de- 
velopment. 

In the realm of reactor power, with the 
exception of the combat submarine systems, 
there exists only one program to develop a 
true deep-submergence nuclear power plant — 
the NR-J project. No other reactors for under- 
sea use have been developed to the proof-of- 
principle stage. 

In addition to the need for a definition of 
oceanographic objectives, the basic obstacle to 
nuclear ocean power is cost. Less familiar 
is the problem of designing nuclear hardware 
to operate in the marine environment. The 
solution to both depends on resolving the old 
problems of system efficiency, compactness, 
and weight—especially shield weight. As the 
need for nuclear ocean systems becomes more 
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evident, the incentive will lead to solutions to 
these problems. 

It is clear that nuclear-powered devices have 
a definite and often undisputed place in the 
exploration and exploitation of the oceans. 
Wherever long-endurance power is required, 
nuclear systems are serious competitors with 
fossil fuel and electrochemical energy sources. 
Although the development of nuclear power for 
the sea is still in a very early stage, the re- 
sults thus far have been so promising, and the 
variety of potential applications of nuclear 
power so great, that continued and expanded 
research and development in this area appear 
nearly inevitable for the future. 
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Advances in Power-Reactor Shield Analysis 


By Phillip S. Mittelman, United Nuclear Corporation* 


Power-reactor shields are massive structures 
designed to reduce external radiation heating, 
radiation dose, and radiation activation to ac- 
ceptable limits. Because the reactor shield 
represents a significant cost and also influences 
the design of the reactor building and contain- 
ment structure, it is worthwhile to optimize 
the shield design so as to minimize both direct 
and indirect costs. 

Toward this end there has been a slow but 
steady advance in the state of the art of power- 
reactor shield analysis. Accurate prediction of 
shield performance is the first step in an 
optimization program—and in itself can re- 
duce costs by eliminating the need for extra 
shield thickness to compensate for calculational 
uncertainties. 

Advances in Shield analysis have consisted 
mainly of transforming feasible methods of 
shield analysis into practical methods by the 
use of computers. Most of the fruitful activity 
has been associated with: (1) improved tech- 
niques for analyzing bulk shields and (2) analy- 
sis of complex geometries by means of the 
point-kernel or Monte Carlo method. 


Bulk-Shield Analysis 


A substantial portion of a power-reactor 
Shield can be regarded as a one-dimensional 
(slab, sphere, or cylinder) multilayer assembly 
of material surrounding a source region of 
similar geometry. A shield-analysis technique 
attempts to determine the neutron and gamma- 
ray leakage from such a shield, taking account 
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of secondary gammas produced by neutron 
interactions in the shield. In addition, the 
radiation heating and activation of shield com- 
ponents are also desired outputs. 

Several very useful approaches have been 
developed in recent years to study bulk shields 
quickly and accurately. In general, these meth- 
ods have replaced the older techniques that 
involved the performance of power and geo- 
metrical transformations! on experimental bulk- 
shield configurations. They are more accurate 
and can be performed for situations in which 
experimental data are not available. Because a 
complete bulk-shield analysis can be completed 
in less than 15 min of IBM-7090 time, the 
effects of design modifications can be studied 
easily. 

The most widely used methods combine a 
removal theory and multigroup diffusion theory. 
The United Kingdom RASH codes,” the Swedish 
NRN method,’ and the Hanford MAC program‘ 
divide the high-energy region into a number of 
energy groups and assume that source neutrons 
are attenuated exponentially. Energy-dependent 
removal cross sections are employed; neutrons 
“removed” from the high-energy groups slow 
down in the shield and lead to secondary gamma 
sources. In the RASH procedure, continuous 
slowing down is used, whereas NRN and MAC 
permit slowing down from one group to any 
other group. In these programs the capture- 
gamma sources are determined, and then 
gamma-ray transmission is calculated to de- 
termine the leakage gamma flux. 

The RASH and MAC programs calculate 
gammas by point-kernel techniques with ap- 
proximate multilayer buildup factors. However, 
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NRN uses a highly efficient slab-geometry 
Monte Carlo program that provides gamma 
heating and leakage data. Figures 1 to 3 show 
some typical results obtained with the RASH, 
MAC, and NRN codes. When these methods are 
applied to water-reactor shields consisting of 
water, iron, and concrete, they can predict 
fast-neutron and gamma doses that are accurate 
within a factor of about 3. Both RASH and NRN 
have been coded for Ferranti mercury com- 
puters; MAC is a FORTRAN program that is 
usable on the IBM-7090—7094 series and canbe 
adapted easily to other large computers. The 
latest version> of NRN is programmed in 
FORTRAN IV and is operational on the IBM- 
7044. The latest version of RASH, RASH-F, 
is being programmed in FORTRAN and should 
soon be available. 


MULTIGROUP TECHNIQUES 


It is tempting to use the multigroup diffusion 
and transport programs* ® familiar to reactor 
analysts to perform the shield analysis. Such 
programs are useful to determine capture- 






























10 2PPvg 
ay WOT KH HHH ne 
1912 10 | 
109° ee: | 
800 900 
1910 
- Air gap 
& 108-- 
+ Air gap Hetaestamite 
2 106 Pressure 
= vessel Concrete 
= ihe f | (biological shield) 
3 Thermal Air gap 
104 shield 
S 
2 | 
“a AUS 
2 s 
= 
© , RASH-E 
109}— \v 
10“ 
b—— 
AP. ey eee eee ems ee 
500 600 700 800 900 
Distance from Core Axis, cm 
Fig. 1 Thermal-neutron flux distribution through 


yadial bulk shield of No. 1 reactor at Calder Hall op- 
erating at 180 Muwu(t) showing values calculated by 
RASH codes and measured data points.?” 
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Fig. 2. Gamma distribution through Hanford Shield 
Facility showing points measured® and curves cal- 
culated by MAC code.4 


gamma sources in the inner shield layers but 
are inadequate to handle deeper penetration. 
Generally, high-energy radiation scattered 
through small angles is responsible for the 
dose at the outer shield layers, and reactivity 
programs underestimate it. Analysts have ob- 
tained useful shielding results by performing 
multigroup calculations with P, theory or the 
Carlson S, methods and then normalizing the 
high-energy fluxes or doses to the results of 
accurate high-energy calculations, i.e., mo- 
ments method or Monte Carlo method. The 
P,MG method’ uses this technique, and Tomonto”” 
has reported on a similar calculation combining 
the Monte Carlo and S, methods. 

The P;MG program" uses a P; rather than a 
P, approximation to the transport equation and 
appears to be accurate at deep penetration 
(~150 cm) in slab shields without needing 
normalization to more refined methods. 

Special versions’”*® of the Carlson S, pro- 
gram have been prepared for shielding applica- 
tions. These programs can incorporate highly 
anisotropic angular distributions following elas- 
tic scattering and thus are able to treat more 
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Fig. 3 Comparison of NRN, NIOBE, and RASH-B calculations with measurements for two shield configu- 
rations.” Note that all neutron fluxes are presented as ratios of those data to the NRN prediction. Gamma 
data are presented as ratios of exposure rates to the RASH-B prediction. 


accurately the high-energy neutron-transport 
problems. Although they have not yet been 
applied to full power-reactor shields, they have 
been used successfully for limited depths of 
penetration, and DTF-IV has been successful in 
calculating gamma-ray penetration. 

The NIOBE program (Numerical Integration 
of the Boltzmann Equation)'* represents per- 
haps the most ambitious attempt at the numeri- 
cal solution of the Boltzmann transport equation 
for neutrons and gamma rays in spherical 
multilayer assemblies. It has been used suc- 
cessfully in Europe to study power-reactor 
shields and to check out the RASH’ and NRN*® 
methods. The requirements of 2.5 hr of IBM- 
7094 time and the absence of a thermal group 
have restricted its application for water-reactor 
shields. 


Complex-Geometry Problems 


In power-reactor shielding the designer often 
is faced with the problem of determining the 
radiation leakage or streaming through shield 


penetrations, such as coolant lines or ducts, 
instrument tubes, and gaps around shield plugs. 
Shield heating problems often require that the 
energy deposition in particular locations in the 
inner shield be determined fairly accurately 
(within 20 to 30°c). In these instances a one- 
dimensional representation of the shield is not 
adequate, and two- or three-dimensional repre- 
sentations are needed. Generally such problems 
are handled by the use of one of the following 
techniques: (1) point kernel, (2) combination 
point-kernel—multigroup, or (3) Monte Carlo. 


POINT-KERNEL METHODS 

The point-kernel programs break the source 
region into a number of equivalent point sources 
and determine the contribution of all such points 
to one or more detectors. The contribution of 
any point to any detector is a function of the 
distance between source and detector and the 
amount and type of material on the line joining 
them. A “kernel” predetermined from experi- 
ment or more exact calculation provides the 
material attenuation information. 
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The most useful application of point-kernel 
programs is for estimating fast-neutron dose 
or gamma dose. For fast-neutron dose the 
method is primarily used for shields containing 
a thick final layer of hydrogenous material. 
The attenuation of other layers is accounted 
for by the use of removal cross sections,'® and 
a special kernel is used for the hydrogenous 
layer.'’ Nonhydrogenous shields (carbon and 
beryllium) have been analyzed with some suc- 
cess by the use of kernels derived from mo- 
ments-method calculations.'® 

For gamma-ray dose calculations, the kernel 
is the product of the uncollided dose and a 
buildup factor. The latter is taken to be the 
buildup factor appropriate for the last layer, 
or one of several multilayer-buildup-factor 
formulas is used.'*-?° Typical programs are the 
GE-14-0 program’! and the Los Alamos QAD 
program.” Each of these programs handles 
complicated figures of revolution. 

When the radiation actually travels along the 
rays in question, point-kernel methods can give 
accurate answers. At times, however, the de- 
tector response is predominately caused by 
radiation that has scattered around local shield- 
ing. Then the point-kernel methods can be 
orders of magnitude wrong. Each problem 
must be examined to see if such possibilities 
exist. If so, then Monte Carlo techniques should 
be used. In fact, it would be good practice to 
check a point-kernel calculation for a particular 
shield against Monte Carlo or to experiment 
before this technique is used extensively. 


COMBINATION POINT-KERNEL—MULTIGROUP 
METHODS 

The various kernels discussed above gen- 
erally yield information on dose, aeating, or 
some other integral quantity. {f more detailed 
spectral information is requir-d, it generally 
can be obtained when (1) a one-dimensional 
problem is solved with an approximate repre- 
sentation of the shield layers present between 
the core and the detector in question and (2) 
the fluxes are scaled so that the computed dose 
equals the dose determined from a point-kernel 
calculation for the same detector point. The 
RASH, NRN, MAC, P,MG, P3MG, or DTF cal- 
culations can be used for the one-dimensional 
analysis. 


MONTE CARLO METHODS 


When accurate answers for complex geom- 
etries are required, the Monte Carlo programs 
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can be used. These programs track radiation 
through space and simulate its actual interac- 
tions with matter. To the extent that good cross 
sections are available, these programs can 
give highly accurate answers to the shielding 
problem. 

Two factors must be kept in mind when a 
complex-geometry Monte Carlo program is 
used: (1) if improperly used, the program will 
give answers withan apparently small statistical 
uncertainty but with a large absolute error, and 
(2) the running time for problems may be of the 
order of hours on machines in the IBM-7094 
class. 

The potential inaccuracy arises from the 
possibility that important spatial and energy 
radiation sources are not properly sampled. 
For instance, the neutron dose outside a thick 
water shield is produced largely by those neu- 
trons whose initial energies are greater than 
6 Mev. The Monte Carlo program may pick its 
source neutrons in the reactor core from a 
fission-spectrum energy distribution. If 1000 
particles are being followed, the probability 
of picking a neutron above 6 Mev is So low that 
only approximately 20 will be picked. If these 
20 are picked isotropically, few, if any, may 
actually enter the shield. In this situation a 
low emergent dose will be calculated, and the 
code user will be unaware of the error. An 
experienced shield physicist or engineer would 
certainly run a special high-energy-source 
problem to ensure that the important source 
energies are sampled properly. 

There are several complex-geometry neutron 
and gamma-ray Monte Carlo programs in routine 
use today. The principal programs are: 05R,”° 
FMC,”4 ADONIS,”® and UNC—SAM.”° 

Some of the important characteristics of these 
programs are tabulated below. Because both 
O5R and UNC-—SAM still are growing in capa- 
bilities, it is suggested that potential users 
contact the originating group to determine the 
applicability of the program to their specific 
problem. 

Any organization contemplating the use ofa 
complex-geometry Monte Carlo program should 
allow a learning period of several man-years 
before the program will be operational. These 
programs should be regarded more as experi- 
mental shielding facilities than as computer 
programs. They require skilled scientists or 
engineers for operation and continuous main- 
tenance and adjustment by programmers to fit 
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new problems. These programs should be kept 
up-to-date as new techniques of importance 
sampling and routines for representing nuclear 
interactions are developed. 


O5R 


Computer: Operational on CDC-1604A and IBM-7094, 

Language: Combination of FORTRAN and machine 
language. 

Geometry capability: Handles a three-dimensional 
array of figures bounded by planes and quadric 
surfaces. 

Radiations: Neutrons. 

Output: Produces interaction tapes that can be edited 
by user to evaluate fluxes, secondary gamma 
sources, heating, etc. 

Variance reduction methods: Splitting or Russian 
roulette on region weights as functions of energy 
and position. 


UNC-—SAM | 


Computer: Operational on CDC-1604A and Aberdeen 
BRLESC. 

Language: FORTRAN. 

Geometry capability: Handles a three-dimensional 
array of bounding boxes, each of which can contain 
up to seven boxes, spheres, cylinders, or wedges. 
Each of the latter can contain up to seven figures, 
which in turn can contain seven figures and so on. 
Orientation of figures inside bounding boxes is 
arbitrary. 

Radiations: Neutrons and gamma rays. 

Output: Neutron and gamma fluxes in each region or 
at specified points as functions of energy. Interac- 
tion tape yields secondary neutron and gamma 
sources and serves as source tape for secondary 
problem. Response functions are dose, heating, 
etc. Transmitted particle tape contains phase space 
coordinates of radiation leaving emergent face of 
configuration and can be used to continue problem 
into additional regions. 

Variance reduction methods: Splitting or Russian 
roulette on region weights as functions of energy 
and position. 

Cross sections for neutrons: 18 Mev to 0.037 ev (80 
enc: gy intervals) plus one thermal group; isotropic 
or anisotropic elastic scattering; level or con- 
tinuum inelastic scattering; absorption; (”,2n) re- 
actions in beryllium. Gamma rays attenuated by 
Compton scattering and absorption by photoelec- 
tric and pair-production processes. 

Special features: Subroutines for preparation of 
cross-section data from point data and to simplify 
geometry input. 


ADONIS 


Computer: Operational on IBM-7090-— 7094. 

Language: FORTRAN and FAP. 

Geometry capability: Handles a three-dimensional 
array of boxes. 

Radiations: Neutrons and gamma rays. 

Output: Neutron and gamma fluxes in each region as 
functions of energy. Interaction tape yields secon- 
dary gamma and neutron sources and serves as 
source tape for secondary problem. Response func- 
tions are dose, heating, etc. Transmitted particle 
tape contains phase space coordinates of radiation 
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leaving emergent face of configuration and can be 
used to continue problem into additional regions. 
Variance reduction methods: Splitting or Russian 
roulette on region weights as a function of posi- 
tion only. 
Cross sections for neutrons: Same as UNC-—SAM 
except that no thermal group is present. 


FMC-N and FMC-G 


Computer: Operational on IBM-704 and IBM-7090. 

Language: Machine language. 

Geometry capability: Handles a three-dimensional 
array of figures bounded by quadric surfaces. 

Radiations: Neutrons and gamma rays. 

Output: Radiation fluxes in volumes; leakage tallies; 
energy deposition. 

Variance reduction methods: Splitting or Russian 
roulette on region weights as functions of energy 
and position; exponential transformation; statistical 
estimation of fluxes and leakages. 

Special features: Allows reflecting boundaries to 
simplify geometry description. 


Conclusions 


Techniques are available to study water- 
reactor bulk shields rather quickly and with 
accuracies of factors of about 3 for the neutron- 
and gamma-leakage doses. These techniques 
are embodied in a variety of computer pro- 
grams which, with the exception of P;MG, use 
a transport kernel for the high-energy neutrons 
and a continuous slowing-down or multigroup 
slowing-down model to determine low-energy 
fluxes and capture gamma sources. 

If geometrical considerations are important, 
approximate answers can be obtained by theuse 
of point-kernel programs. These operate in a 
few minutes of IBM-7090 time but are of un- 
known accuracy. They should be normalized to 
either experimental results or Monte Carlo 
results in the same geometry. 

Monte Carlo methods are available for the 
complex-geometry problems associated with 
ducts, plugs, and structural heating. When used 
properly, these programs are capable of ac- 
curacies of factors of 2 for neutron and gamma 
doses. 
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Measuring Neutron-Flux Distribution 
in Operating Reactors 


By Glenn F. Popper 


Reactor economics, safety, or experimental 
control often provide the incentive for accu- 
rately and continuously measuring the distribu- 
tion of neutron flux within a reactor core. Thus 
designers of both power and research reactors 
find it worthwhile to consider instrument sys- 
tems for making such measurements. ! 


Reasons for Monitoring Flux 


Greater power, burnup uniformity, safety, 
and radiation monitoring are among the pur- 
poses for measuring flux distribution. 

Core uprating is an economic reason for 
monitoring the flux in a power reactor. The 
flux-mapping system used in the Yankee Rowe 
reactor, for example, confirmed a conservatism 
of at least 20% in the calculated maximum-to- 
average heat fluxes under the most unfavorable 
control-rod conditions.? On the basis of this 
information, the AEC authorized Yankee to 
increase plant gross output from 392 to 485 
Mwit) in June 1961. Power operation of core 2 
began Sept. 21, 1962. Again, as happened in 
core 1, the in-core instruments confirmed that 
higher powers could be achieved. Throughout a 
significant portion of core 2 life, the reactor®»4 
operated at 540 Mw(t). Core 4 is expected 
to be able to operate at 600 Mw/(t) solely 
because of the information obtained from the 
flux-monitoring system. 

The economic success of the Yankee power- 
reactor system stimulated designers to incor- 
porate flux-monitoring systems, of one kind or 
another, in the design of other reactors, such 
as San Onofre, Saxton, and Connecticut Yankee.° 


The Connecticut Yankee reactor, to be opera- 
tional in 1967, will incorporate insertion drives 
similar to the Yankee flux-monitoring system 
but will substitute miniature ion chambers for 
activation wires.® These detectors were satis- 
factorily tested in the Saxton reactor in 1964. 
However, the Saxton reactor routinely uses a 
flux-wire system similar to the Yankee Rowe 
system. San Onofre is scheduled to use an 
“Aeroball” system’ (described later). 


Another example of the successful use of a 
flux-monitoring system to economic advantage 
was in the Dresden 1 reactor.® As a result of 
work with in-core ion chambers, the Dresden 
core was actually taken to rated power [630 
Mw(t)} with 7.5% of the fuel removed. These 
examples show that the ability to determine 
experimentally, and with reasonable accuracy 
(~5%), the local flux distribution can create 
long-term economic savings that more than 
compensate for the high initial cost (usually 
greater than $50,000) of the instrument system. 
Other reactors that are using in-core flux- 
monitoring systems are Big Rock Point, Hum- 
boldt Bay, SENN, and the Japanese Power 
Demonstration Reactor.® 


Uniform fuel burnup also can lead to sub- 
stantial savings in either power or research 
reactors. With the information derived directly 
from the flux-monitoring-instrument system, 
the reactor operator can usually adjust the 
control rods to flatten the shape of the flux 
across the whole core and thereby obtain a 
more uniform burnup that will result in a 
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smaller fuel inventory and a lower fuel- 
management cost. 

From the standpoint of reactor design, the 
engineer or physicist can better evaluate the 
core performance from fuel and thermal con- 
siderations if flux maps are available.'” For 
example, wire irradiations in the SENN core 
have provided data that :can be used either 
to verify present calculational methods or to 
provide correlation factors for use in future 
reactor designs. In this way, advanced designs 
can be realized without extreme conservatism. 


SAFETY CONSIDERATIONS 


Safety in large power reactors can also be an 
important reason for requiring in-core flux 
monitoring."! The Dresden nuclear power plant 
contains an operational in-core monitoring 
system of 64 points as required by the AEC 
license. This system is connected directly to 
the safety logic circuits. All future large 
prototype power reactors will probably be 
required to have a flux-mapping system con- 
nected to the safety system. 

In the Halden Boiling Water Reactor (HBWR), 
core power was limited to 20 Mw/(t) because 
of expected hot-channel factors. In-core flux 
measurements made with gamma thermom- 
eters verified that, with only a few in-core flux 
monitors connected to the safety circuit, an 
increase in power to 30 Mw(t) could easily be 
justified.’ This 50% power increase, accord- 
ing to the report, represents almost $4 million 
per year in additional revenue. 


IRRADIATION MONITOR 


The experimenter using a test or research 
reactor for irradiation studies is almost al- 
ways faced with two serious problems. First, 
because of the very frequent changes in the 
experimental equipment installed within the 
reactor and the subsequent control-rod-position 
changes required or because fuel-loading 
changes are necessary, an intermediate mea- 
surement of the neutron flux available at any 
particular irradiation location is not obtain- 
able. The second problem is that, even after the 
flux has been determined, continuous changes 
occur with time. Hence installation of a flux- 
monitoring system is one solution to these 
problems. As a matter of convenience and 
economy, it would be preferable for the reac- 
tor facility to provide such a system for 
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continuous neutron-flux-distribution measure- 
ments. However, no research or test reactor 
that will provide such continuous measure- 
ments is known. 


Measurement Techniques in Use 


Now that the reasons for making these mea- 
surements have been considered, it is appro- 
priate to review some of the techniques with 
which to make neutron-flux-distribution mea- 
surements. By far the most familiar and 
popular method involves activation techniques. ! 
This method uses any material that can detect 
neutrons by being measured for induced radio- 
activity after exposure in a neutron environ- 
ment and is classed as a steady-state technique 
giving integral flux measurements. Materials 
that have well-defined absorption cross sec- 
tions are usually used. Gold, indium, manga- 
nese, copper, cobalt, sulfur, and phosphorus 
are some common solid materials used. How- 
ever, liquids or gases could also be used. The 
solid materials take the shape of foils, wires, 
pellets, or spheres; most are available com- 
mercially. '4 


FOILS 

Counting techniques are readily available 
for foil-activity determinations, and foils are 
being successfully used in many reactors for 
neutron-flux measurements. For example, 
gold and indium foils have been installed in the 
Oak Ridge Graphite Reactor for neutron-flux 
determinations.'* Similar measurements" using 
both foils and wires were made in EBR-II. The 
activants used were 7*U and 7°y, 

Wires and spheres or pellets are gaining 
favor over foils as a solid activant shape. Two 
of the reasons behind this trend are that they 
can be inserted into and removed from a reac- 
tor core more easily, even during reactor 
operation, and the commercial availability of 
complete wire- and ball-drive systems (Refs. 
2, 7, 9, 18). 


WIRES 

A wire technique is used successfully in the 
Yankee reactor.? Activation wires made of 
carbon steel with approximately 0.5% manga- 
nese can be installed in any of 22 core posi- 
tions, as shown in Fig. 1. Three separate 
irradiations of the eight available flux wires 
are required for a complete core mapping. 
After the wires have been irradiated about 
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Fig. 1 Yankee in-core instrumentation showing 22 
core positions into which activation wires can be in- 
serted.? 


4 min at 392 Mw(t), which corresponds to a 
time-integrated thermal flux of approximately 
4.3 x10" neutrons/em’, they are withdrawn 
from the core and stored approximately 1 hr 
while the short-lived radioactive isotopes de- 
cay. The gamma activity from Mn is then 
scanned with a Nal scintillation crystal— 
photomultiplier-tube combination. The tube 
current output is displayed directly on an x—y 
plotter (x is axial core position; y is neutron 
flux), but it is also fed to an integrator where 
the integrated flux is plotted against axial core 
position on another x—y recorder. 

Activity buildup in the wires in numerous 
irradiations, **Mn decay during the elapsed 
time between irradiation and counting, wire 
inhomogeneity, and activation during insertion 
and removal require that corrections be ap- 
plied to the data during reduction. These nec- 
essary corrections complicate the analysis and 
make fully automated data-reduction equipment 
expensive. 

Complicated and often unique mechanisms 
are required for traveling-wire activation and 
analysis.'? The typical mechanism shown in 
Fig. 2 was used in the Vulcain core experi- 
ments in the BR-3 reactor. Such a device adds 
Significantly to the cost of the entire flux- 
monitoring system. The accuracy to which a 
wire is positioned is limited by the flexibility 
of the manipulation cable in the thimble and can 
vary by as much as 1 in. Silver or stainless- 
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Fig. 2 Typical mechanism for traveling-wire acti- 
vation and analysis showing general arrangement of 
flux-plotting machine” for the BR-3 reactor. 


steel neutron absorbers have been fixed to the 
detector thimbles at the top and bottom of the 
core to overcome these positioning errors 
(by providing fixed reference points on the 
wire). Position accuracy,” determined by the 
absorber-placement accuracy, is 0.06 in. 


SPHERES 


A spherical shape is also being used for 
activation and has certain advantages over 
wires."’!?»21 The very small sphere has suffi- 
ciently low radioactivant concentration to keep 
self-shielding low and approaches the ideal 
point detector for isotropic flux determination. 
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This shape eliminates obliquity variations and 
results in insignificant flux depression and 
perturbation factors. Because the spheres are 
inserted into and withdrawn from the reactor 
core pneumatically, there is no need for the 
expensive and bulky shielding that is necessary 
for the highly activated portions of the manipu- 
lation cable used with wire-drive systems. The 
pneumatic system, however, must be vented 
to an exhaust stack. Because of the positive 
thimble stops, the chain of spheres can be po- 
sitioned very inexpensively, quickly, and more 
accurately in the core than can wires. The 
spheres are also particularly adaptable to a 
cheaper automatic-analysis technique because 
the system can count individual spheres. Cor- 
relation between each sphere and core position 
is well defined and limited only by the mea- 
surement accuracy of the ball or head stops. 
However, when individual spheres are counted, 
a direct indication of relative flux distribution 
is not available. 


EBWR EXPERIENCE 

Small, approximately 0.063-in.-diameter, 
aluminum—1% cobalt spheres (Fig. 3) were 
inserted pneumatically into six tubes located in 
selected core positions in the EBWR.”! The 
spheres were encapsulated in small pieces of 
stainless-steel tubing, 0.096 in. in inside diam- 
eter, by forming the tubing over the end beads 
of 0.058-in. commercial stainless-steel bead 
chain, and thereby provided an easy means for 
accurately spacing the activant spheres. One 
disadvantage was that the capsules had to be 
machined open with a small lathe in a glove 
box to retrieve the spheres. This problem 
could be overcome by the use of aluminum or 
plastic tubing that is opened with a simple 
cutter. A sample changer was used so that in- 
dividual spheres could be pneumatically trans- 
ferred to a detector assembly, counted, and 
then returned to an indexed storage location in 
the changer. A multichannel analyzer with 
automatic printout was used for routine analy- 
Sis of each sphere. Plans to use a slightly 
modified version of the system for the Pluto- 
nium Recycle Test Program are being made. 


AEROBALL SYSTEM 


More recently Westinghouse has developed 
the “Aeroball” system, in which in-core flux 
monitoring is accomplished by the pneumatic 
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Fig. 3 Stainless-steel tubing encapsulates alumi- 
num—1% cobalt spheres inserted into six tubes 
located in selected EBWR core positions.*! 


movement of a column of separate 0.063-in.- 
diameter balls into the core for irradiation and 
then to a scintillation detector for scanning. 
This system is being used in the SELNI reactor 
in Italy and in the San Onofre reactor."? The 
Aeroball thimbles, shown in Fig. 4, enter the 
San Onofre vessel through five spare control- 
rod ports and are then routed into selected 
core positions. The ball set used in each 
thimble is irradiated, put into storage, and 
then moved individually by diverter valves to a 
rotary readout table for detection. The ball 
set to be scanned enters the table and is 
formed into a horizontal circle around the 
counting table. The table then rotates past the 
scintillation detector and the 1-in.-wide colli- 
mator; 16 balls are scanned at a time, and the 
integral flux is plotted. Results from the SELNI 
system show that a more detailed flux profile 
can be obtained with a narrower collimator. 


GAS ACTIVANTS 


Gases are also being used as activants.” 
Argon is an excellent monitor of thermal- 
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Fig. 4 Aeroball system for monitoring flux uses 
0.063-in.-diameter balls that are moved pneumati- 
cally in thimbles within core, irradiated, then scanned 
by scintillation detector.’ This system will be used 
in the SELNI and San Onofre reactors. 
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neutron flux. The only significant isotope, 
“Sar (99.6%), has a 0.53-barn thermal-neutron- 
absorption cross section with no resonance 
peaks. The other isotopes do not decay with 
gamma emission after neutron capture. Neu- 
tron capture transmutes ‘Ar into “‘Ar, which 
decays by emitting a 1.29-Mev gamma photon 
in 99% of the disintegrations. Since argon de- 
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cays to stable daughters, it is not a source of 
contamination and can also be used as an inert 
protective atmosphere. 

This argon-activation method was used to 
measure neutron flux in an irradiation experi- 
ment in the ORR. The argon gas flows in and 
out of the irradiation capsule through con- 
centric tubes that are small compared to the 
capsule. This size difference produces a much 
higher gas flow rate through the tubes, com- 
pared to the flow through the capsule, and 
causes a lower activation in the tubes for a 
given flux density. Since most of the tubes are 
located in a much lower neutron flux than the 
capsule, all the activation can be assumed to 
occur within the capsule. However, since ex- 
posure time is sensitive to capsule tempera- 
ture, an empirical correction must be made. 
The delay time between activation and detec- 
tion, usually only a few seconds, depends on 
the gas flow rate and the flow path length 
between capsule and detector. The largest 
error in flux determination is attributed to the 
approximate 10% accuracy of the gamma- 
spectrometer detector. Moving the capsule in 
the irradiation thimble permitted neutron-flux 
distributions to be obtained. 


NEUTRON SPECTRUM 


The various analytical methods used to de- 
termine a neutron spectrum from activation 
measurements range from simple calculations, 
requiring minimum cross-section information, 
to complex computer techniques, requiring de- 
tailed information that reproduce the spectrum 
with high-energy resolution. The five basic 
methods currently being considered are:'* (1) 
a resonance-reaction method, (2) an age-mea- 
surement method, (3) an effective threshold- 
reaction method, (4) a polygonal or polynomial 
method, and (5) an orthonormal-expansion 
method. In all cases the selection of a specific 
radioactivant material to define a desired neu- 
tron energy is a great problem. Additional 
work is needed in this area of material selec- 
tion, especially for measurements to be madé 
in the 10°-ev neutron region. 


MINIATURE ION CHAMBERS 


Although activation techniques will continue 
to be necessary for neutron-spectrum mea- 
surements, miniature ion chambers are rapidly 
replacing activants as in-core flux-distribution 
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detectors. These ion chambers are small 
(some with diameters <0.090 in. and only 1 in. 
long) and have linear response and continuous 
outputs. Nearly the entire flux range [from 
about 10‘ neutrons/(cm?)(sec) up to about 10'4 
neutrons/(cm’) (sec)| of today’s cores can be 
monitored by in-core ion chambers.”° 


Fixed-position in-core ion chambers, as 
Shown in Fig. 5, have been used to monitor 
neutron-flux distribution in the Dresden 1 
reactor (Refs. 8, 11, 20, 24, 25). The Dresden 
in-core power-monitoring array consists of 16 
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Fig. 5 Fixed-position in-core ion chambers of the 
type used to monitor neutron-flux distribution in the 
Dresden 1 reactor.” 
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assemblies for determination of radial dis- 
tribution, each of which includes four ion 
chambers distributed at discrete elevations to 
determine axial distribution. A dry thimble to 
accept an activation wire for calibration pur- 
poses is also included in each assembly. The 
chambers are calibrated at low power to cor- 
rect for variations in the “as-built” sensitivity. 
Calibration is accomplished by the insertion 
of a copper—titanium wire in each of the 16 
ion-chamber guide thimbles during operation 
at approximately 50 Mw(t) and by irradiation 
for about 3 hr. The relative calibration of each 
ion chamber is then attained when the in- 
dividual outputs are corrected to give a flux 
distribution equal to that obtained from the 
irradiated-wire distributions. Thereafter ac- 
curacy of about 5% is achieved. Calibrations 
are made as needed for correction of burnup of 
the fissionable uranium coating used in the 
chamber to render it neutron sensitive. The 
coating is an aluminum—20 wt.% fully enriched 
uranium alloy. The life of such a chamber is 
relatively short; 50% uranium-coating burnup 
corresponds to an integrated neutron flux of 
approximately 107! neutrons/ cm’, approximately 
nine months in the Dresden reactor.”° 

Traversing ion chambers are rapidly be- 
coming the most popular means of monitoring 
neutron-flux distribution in reactor cores.’ 
Much longer lives are achieved because of the 
shorter core-exposure time. A system that 
uses activation wires for low-power operation 
[~10° to 10! neutrons/(cm’)(sec)] and travers- 
ing ion chambers for high-power operation 
[~10" to 10'* neutrons/(cm’)(sec)] was planned 
for the Experimental Gas-Cooled Reactor.”4:?6 
The ion chambers were to be used as the pri- 
mary means of monitoring neutron flux in the 
core. The chamber-cable probes would be in- 
serted into 18 thimbles, and a flexible multiple- 
strand cable would serve to position the cham- 
bers in the core. The center of this cable would 
be filled with a triaxial electrical lead. The 
outer wire strand of the cable would be an open 
helix and would mesh with a geared drive wheel 
either manually or electrically driven. 

The chamber components and wires of the 
cable would be made of commercially pure 
titanium; quartz fiber was specified as elec- 
trical insulation in the cable. Forsterite would 
act as an electrical insulator in the chamber 
which was lined with '°B and would have a 
thermal capability of 1200°F. This ‘°B coating 
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was selected to minimize the radioactivity of 
the chamber assembly after withdrawal from 
the core. The relatively rapid burnup of the 
boron (33% at 10° neutrons/cm”) is acceptable 
because the probes would be used for periodic 
traverses only. Recorders in the control room 
plot neutron flux against position along the path 
of travel through the core. The positioning 
accuracy would be the same as with wire-drive 
systems. 


THERMOMETERS 

Gamma and neutron thermometers have also 
been applied to the measurement of neutron 
flux. 19.20 By this method, measurement is made 
of the temperature rise above ambient caused 
by fission or gamma heating. Several practi- 
cal devices have been developed, but none 
are known to be available commercially.?’-% 
Nevertheless, thermometer type sensors con- 
tinue to be alluring as detectors. The obvious 
advantages of basic simplicity, easy construc- 
tion, and simple readout systems are joined 
with the possibility of making a detector that 
physically resembles the fuel element. This 
could result in a power monitor, in contrast 
to a neutron-flux monitor. The problem of 
analytically relating the neutron flux to local 
power could then be eliminated, if local power 
actually is the desired parameter. 

The neutron thermometer’’ used by Du Pont 
at Savannah River is shown in Fig. 6. This 
consists of a thermocouple with a uranium— 
aluminum button welded to the hot-junction end. 
The heat generated in the button increases the 
temperature of the sample thermocouple above 
that of the reference. The difference signal is 
proportional to the neutron flux. Although this 
thermometer has shown a sensitivity of about 
5°C per 10” neutrons/(cm?)(sec), it is limited 
by the uranium burnup. Sensitivity would be 
reduced about 10% at 5x10”? neutrons/cm’. 
The time response is about 0.5 sec. 


SOLID-STATE DEVICES 


Solid-state devices*!;** have proved useful 
as compact neutron detectors for measuring 
flux distribution in critical experiments and in 
reactors in which the flux is in the range of 
10‘ to 10! neutrons/(cm?)(sec). Operation in 
high flux is impractical because a dose of 
10? neutrons/em? damages the semiconductor 
crystal lattice. For example, the device has a 
life of only 100 sec in a flux of 10" neutrons/ 
(cm?)(sec). The operators overcame this diffi- 
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Fig.6 Neutronand gamma thermometers. Heat gen- 
erated in uranium—aluminum button raises tempera- 
ture above that of reference thermocouple, and the 
neutron flux is determined from the difference signal, 
which is proportional to the neutron flux.2? Gamma 
thermometer measures temperature rise above am- 
bient caused by gamma heating. 


culty by removing the semiconductor from the 
reactor core but leaving a fissionable target 
within an evacuated tube that extends into the 
core region.*® The detector then counts fission 
fragments from the target, and the neutron-flux 
distribution can be plotted if the target and tube 
are moved within the core. This method mea- 
sures neutron flux directly within 3%, and 
the target position is accurately known. This 
method can be used in fluxes up to approxi- 
mately 10'° neutrons/(cm’) (sec). 


Advanced Techniques 


All the previously discussed devices or sys- 
tems have been or are now being used. 

Self -powered neutron detectors are being de- 
veloped for advanced systems.***’ The neutron 
detector is made in the form: of a coaxial cable, 
as shown in Fig. 7. The current produced in the 
detector comes from neutron-induced radio- 
activity. Detectors using 7"al, ‘In, ‘Ag, and 
other elements have been built, and operational 
data are available. With a rhodium emitter this 
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Fig. 7 Self-powered neutron detector is constructed 
in coaxial cable form;current is produced when emit- 
ter, in this case vanadium, becomes radioactive by 
neutron capture.” 


detector is predicted to have a life 4 times 
longer than **y-coated detectors and 30 times 
longer than those coated with '°B. However, it 
has only recently been introduced as a com- 
mercially available device.* 


Noise systems. The neutron-induced current 
in an ion chamber is caused by a large number 
of individual pulses. The current is composed 
of the steady-state and the varying or noise 
components. The average noise power, or 
steady-state signal, is proportional to the 
neutron flux, and the average deviation, or 
noise component, is proportional to the square 
root of flux.3"38 A practical system using the 
noise component has been developed.”*® This 
system improves gamma discrimination and 
eliminates d-c cable leakage interference cur- 
rents. Hence it can be used in higher tempera- 
ture and higher flux applications. 


Microwave techniques are also being investi- 
gated as a method for neutron-flux measure- 
ments. **-“4 The basic applicable principle is 
that the propagation characteristics of micro- 
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wave energy change as the waves traverse a 
gas that is ionized by neutron bombardment. 
These effects can be feasibly exploited for a 
practical system. The first suggests that the 
attenuation and/or velocity of propagation of 
microwaves traversing a gas plasma is a 
function of the degree of gas ionization. The 
second predicts that there will be cross cou- 
pling of microwave energy between two wave 
guides piercing a common plasma and is like- 
wise dependent on plasma ionization. The third 
effect simply is that a wave guide surrounded 
by an ionized plasma will have induced in it 
microwave noise signals that can be measured. 
Laboratory measurements have confirmed the 
first two effects under simulated reactor con- 
ditions. Research on these detectors is in the 
early stages, but the results are encouraging. 


Regenerative detectors. The useful life of a 
neutron detector is often limited if activant 
or fissionable materials are used. The useful 
lifetime during which neutron detectors can 
provide accurate, reliable in-core measure- 
ments can be substantially prolonged if the 
sensing element incorporates regenerating ma- 
terials such that new fissile atoms are created 
from fertile atoms as the original fissile 
atoms are depleted.**~*’ Detector life can be 
increased by a factor of over 2.5 under certain 
conditions,*® e.g., in ion chambers using a mix- 
ture of both 7y and *%y as the fissionable 
material. The mixture used will depend on the 
type of reactor being instrumented. In addition 
to extending the lifetime of in-core flux moni- 
tors, regenerating detectors will minimize the 
need to calibrate the system as often as for 
fully enriched ?*U detectors. 


Factors Affecting Method Choice 


Before a designer can select any particular 
method for measuring neutron-flux-distribution 
data, usage must be planned. Since neutron flux 
is three dimensional, measurements must usu- 
ally be made in both the axial and radial direc- 
tions. The reactor design greatly influences the 
choice of an instrument. The instrument system 
must be able to withstand the temperature, 
pressure, neutron-flux level, and medium. Cost 
and available manpower are both important 
factors because the development of new sys- 
tems can be very expensive. Table 1 summa- 
rizes the techniques used for neutron-flux- 
distribution measurements and lists some of 
the limitations and significant characteristics. 
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Table 1 CHARACTERISTICS OF NEUTRON-F LUX-DISTRIBUTION MEASUREMENT TECHNIQUES 
Detector System Approximate Neutron- Detector life For 
available available cost of de- flux range, in 10!3 neutrons/ Axial safety 
Technique commercially commercially tector system neutrons/(cm?)(sec) (cm?)(sec) flux resolution circuits Readout 
Activants 
Foils Yes Yes < $1000— 10,000 103 — 1045 <3 months Limited No Integrated 
Wires 10° — 1015 <1 year Continuous No 
Stationary Yes Yes $1000—10,000 Integrated 
Traveling Yes Yes > $50,000 Delayed 
Spheres Yes Yes < $1000-—50,000 103 — 1015 <3 months Limited No Delayed 
Gas Yes No < $1000—10,000 10% 105 <1 day Limited Yes Delayed 
ponnare Continuous 
Stationary Yes > $50,000 Limited Yes 
Traveling Yes > $50,000 Continuous No 
Ion chambers 
Stationary Yes $1000-10,000 1019-1044 <1 year 
Traveling Yes $1000—10,000 101° — 1044 <1 year 
Noise Yes $1000—10,000 10’? — 1044 <1 year 
Regenerative Yes $1000—10,000 1019 — 1044 >1 year 
Gamma, neutron 
thermometers No <1 year 
Solid-state 2 < $1000 104 — 10!! <1 day 
Self-powered Yes < $1000 10° — 104 >1 year 
Microwave No Experimental >1 year 
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Fluid and Thermal 
Technology 


What’s New in Two-Phase Flow 
By Hans K. Fauske 


Two-phase flow is a key factor in both thermal- 
and fast-reactor design and safety. However, 
because flow problems occur in many industries 
and branches of engineering, it is difficult for 
engineers concerned with two-phase flow to 
keep abreast of the progress reported indiverse 
publications. Thus the Symposium on Two-Phase 
Flow that was held at Exeter, England, June 21— 
23, 1965, provided a needed opportunity for 
about 200 engineers of various specialties from 
different countries to communicate concerning 
the latest developments in this area. 

The symposium concentrated heavily on 
“microscopic” aspects of two-phase flow. Al- 
though a few of the papers presented information 
in a form that can be used in reactor design, 
most of them indicated methods and approaches 
for discovering basic information on two-phase 
flow and heat transfer. More specifically, they 
contributed importantly toward the displacement 
of classical two-phase-flow empirical engineer- 
ing by sound physical and theoretical mathe- 
matical approaches that canbe extrapolated with 
much greater assurance to pressures, flows, 
fluids, and geometries other than those tested. 

A summary of the 40 papers presented at the 
Symposium, which was sponsored by the Univer- 
sity of Exeter, will acquaint reactor designers 
with some of the topics discussed, which in- 
cluded flow regimes, wave mechanics, void 
fractions and holdup, two-phase flowcharacter- 
istics, nucleate and film boiling, entrainment, 
flooding, critical flow, and two-phase flow dy- 
namics. 


Flow Regimes 
Detailed descriptions of the flow regimes and 
the regions at which transition takes place from 


35 


one regime to the next must be known to de- 
scribe adequately the heat transfer and pressure 
drop in two-phase systems. 

An experimental study’ of air —water mixtures 
included observations on flow regimes and their 
transitions (Fig. 1). The slug—churn transition 
and the churn—annular transition were de- 
scribed in detail. 
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Fig. 1 Visually observed flow regimes and transi- 
tions of air—water mixtures,! as compared to transi- 
tion predicted by Nicklin’s theory.© 


VERTICAL FLOW 


In an investigation of adiabatic upward flow of 
gas —liquid mixtures, measurements were made 
of droplet entrainment by the gas core and the 
mean film thickness.” Interface aspects were 
observed, and the ratio between liquid film and 
total liquid flow rates were related to volumet- 
ric kinetic energy of the gas (p,Vj,/2), the 
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liquid flow rate, and liquid physical properties. 
The mean film thickness depended on gas 
volumetric kinetic energy, but not on the total 
liquid flow rate. 

Aspects of mass transfer of entrained liquid 
in vertical upflow of air—water were investi- 
gated’ in a 0.375-in.-bore tube with column inlet 
pressures of 25 to 40 psia. Some conclusions 
were: (1) a certain minimum length of tube is 
required before entrainment starts, (2) no ap- 
parent equilibrium is established even after 
long lengths, and (3) for a particular gas flow 
rate, there is a critical liquid flow rate below 
which entrainment does not occur. Order-of- 
magnitude estimates were made of droplet 
mass-transfer coefficients. 

Flow regimes, entrainment, mass inter- 
change, and liquid deposition are critically im- 
portant to the understanding of the burnout 
phenomenon in high-quality mixtures. Because 
two-component mixtures were used in the tests, 
the results cannot be applied directly. However, 
the data are useful in predicting the forms and 
trends to be expected in burnout data. - 


HORIZONTAL FLOW 


An experimental and analytical study* used 
air and mineral oil to determine how the rate 
of atomization (liquid film), equilibrium ratio 
(between liquid flowing as film and liquid flowing 
as mist), and pressure drop in horizontal mist— 
annular flow is affected by fluid properties and 
pipe diameter in the range 51 to140mm. Atomi- 
zation increased with film flow rate and gas 
velocity but decreased with pipe diameter. For 
the equilibrium fraction, 1— 7,, the correlation 
was established that: 


(1 - u,) = 1.247 + 1.04 In [Re? Re# We" 
x (M,/Mg)* (p,/p,)"| 


where Re, equals p,5W%gD/ig, Re, equals 
p,V%_,D/u,, and We equals p,V2. D/c. 


The film thickness was calculated from mea- 
sured pressure drop and film flow rate. 

A test facility has been built at Queen Mary 
College, University of London, to study the an- 
nular flow regime and the mechanism of dry- 
out of the liquid film.’ The high-pressure- 
steam heat source makes operation possible 
under dry-out conditions. The facility, which has 
a wide variety of instrumentation and permits 
visual observation of test-section flow, will 
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enable experiments that should contribute sig- 
nificantly toward the solution of the burnout 
problem in the net boiling regime. 


Two-Phase Wave Mechanics 


The symposium concentrated heavily on 
“microscopic” aspects of two-phase flow, as 
illustrated by the papers dealing with studies of 
the mechanism of wave formation. Although these 
papers do not furnish design criteria directly, 
they are extremely important contributions 
toward the overall understanding of two-phase 
flows. They are prime examples of the kinds of 
experiments and analyses that are needed to 
improve the technology. 

A theoretical analysis that concentrated 
mainly on the regime of roll waves was pre- 
sented.® The method is intended to predict the 
wave Size, speed, spacing, local film thickness, 
and the corresponding two-phase pressure drop. 
A new experimental facility will be used for 
measurements of the variables predicted by the 
analysis; it will test air—water and water —oil 
mixtures. 
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Fig. 2 Effect of countercurrent air velocity on free- 
surface wave frequencies of a liquid film falling down 
(Re, = 80) a vertical surface, as measured at various 
tube positions.’ 
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The effect of a countercurrent stream of air 
on the free-surface profile of a liquid film fall- 
ing down a vertical surface was described." 
A sensitive light-absorption technique was used 
to measure wave properties such as frequency, 
amplitude, and velocity as functions of liquid 
Reynolds number and vertical distance fallen 
by the film. The effects of countercurrent air 
velocities up to 750 cm/sec were examined; 
Fig. 2 shows the effect of air velocity on wave 
frequencies. 

Measurements of turbulence-intensity cor- 
relation coefficients and local shear stress at 
several gas Reynolds numbers and positions 
from a wavy boundary were reported.® The in- 
fluence of the liquid-film wave behavior on the 
gas phase was simulated with a flexible sheet 
that mechanically created standing waves that 
could be oscillated at adjustable frequencies 
(Fig. 3). Spectral analysis of the measurements 
suggested that one mode by which the interface 
acts is to shift the scale of turbulence and that 
for certain conditions the fluid receives energy 
from the interface, resulting in a shear stress 
in the gas phase that is less than would be ex- 
pected if there were no oscillation at all. 
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Fig. 3 Flexible sheet mechanically forms standing 
waves that can be oscillated at various frequencies to 
simulate influence of liquid-film wave behavior on 
gas phase 8 
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INTERFACIAL CONDITIONS 


The following interfacial conditions were ob- 
served in the horizontal flow of air —liquid mix- 
tures: (1) smooth interface, (2) two-dimensional 
waves, (3) three-dimensional waves, (4) roll 
waves, and (5) atomization.? Laws governing 
pressure variations and shear-stress variations 
over a wavy surface were discussed. Detailed 
descriptions of models explaining the physical 
processes responsible for the transitions from 
one type of interfacial structure toanother were 
presented and supported by experimental ob- 
servations. 


Void Fractions and Holdup 


Efforts to determine void volume in two-phase 
systems received strong attention at the sym- 
posium. 


EFFECT OF GEOMETRY 


Because of the lack of adequate data for in- 
dustrial design and scale-up, the liquid volume 
fraction was investigated for vertically upward 
gas—liquid flow in tubes of 0.4 to 12 in. diam- 
eter.'” The liquid holdup was determined by a 
quick-closing-valve technique. Generally ac- 
cepted methods, like the Martinelli—Nelson’! 
and the Hughmark” correlations, were shown to 
predict the liquid volume fraction inaccurately 
over a wide range of air and water velocities 
and column diameters. New correlations were 
presented that included different relations for 
bubble, slug, and froth—annular flow. 

Experiments with rod bundles of annular and 
circular geometries were used to investigate 
liquid volume fraction in gas—liquid adiabatic 
systems.'° In the experiments, which emphasized 
the bubble and dispersed-annular flow regimes, 
argon, water, and ethyl alcohol were used as the 
working fluids at pressures ranging from atmo- 
spheric to 400 psia. The effects of the liquid 
volume fraction on mass flow rate, quality, 
gas density, surface tension, and equivalent 
diameter were examined. Figure 4 compares 
previous correlations with the one proposed 
on the basis of these experiments: 








1 + 


1.358” (2) 
1 + 0.335 - D> \Pc 


where 7” equals 0.9 + 0.050. This empirical cor- 
relation is within 15% of the data obtained. 
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new Cravarolo and Hassid correlation agrees within 
15% with data." 


MEASUREMENT TECHNIQUES 


A dilatation technique was used in an open 
loop to measure total steam volume fractions 
in heated channels.'* The simplicity of this tech- 
nique makes it easier to use in complicated 
geometries than other techniques (X-ray or 
gamma-ray attenuation). Preliminary data of 
total average steam volume fractions in heated 
round ducts were presented. 

The (y,z) reaction was applied to boiling heavy 
water to measure void volume fractions in the 
region of subcooled and low-quality boiling. 
The investigators explored pressures from 10 to 
50 bars, mass velocities from 50 to 1450 kg/(m?) 
(sec), heat fluxes from 30 to 90 watts/cm’, sub- 
cooling from 30 to 0°C, and steam qualities from 
0 to 15%. Pressure, heat flux, and mass velocity 
effects on void fraction were noted in the sub- 
cooled and net boiling regimes. Acomparison of 
this large amount of data with Bowring’s model’® 
.ndicated large discrepancies, particularly for 
lower mass velocities. 

Changes in acoustical properties were dis- 
cussed as another method for measuring voids 
in liquid—gas systems.'’ Expressions were 
derived for the speed of sound and damping in 
bubble mixtures. Figure 5 shows the possibility, 




















ity in aiy—water mixtures at atmospheric pressure 
can determine the vapor volume fraction.!” 


in principle, of determining the volume fraction 
of vapor in mixtures from measurements of 
sound velocity; preliminary experiments with 
air—water mixtures at atmospheric pressure 
prove the technique feasible for very small void 
fractions. 


Two-Phase Flow Characteristics 


Theoretical explanation and satisfactory pre- 
diction of two-phase flow necessitate detailed 
descriptions of the flow processes within the 
system. A number of papers dealt with these 
problems. 

A study of cocurrent vertical mercury —nitro- 
gen flow included measurements of local liquid 
velocities and gas fractions.'® Power-law con- 
tinuity expressions were derived, and velocity 
and void profiles were given for the Hg—N, 
system; Fig. 6 shows typical void-fraction pro- 
files. Levy’s mixing-length theory’? was com- 
pared with computed intensities of density and 
liquid velocity fluctuations. 

An electrical-resistivity probe was used to 
detect the presence of bubbles ina detailed study 
of the statistical properties of two-phase flow, 
primarily in the bubble-flow regime.”° The study 
showed that the order of the Erlangian distri- 
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bution of the bubble-arrival intervals at a given 
point is a statistical criterion for flow-pattern 
identification. 

A theoretical model was presented to de- 
scribe the heat-transfer and fluid-flow charac- 
teristics in vertical boiling channels.”! The equa- 
tions for continuity, momentum, and energy and 
the equation for diffusion and growth of bubbles 
were formulated for the unsteady-state case. 
Radial distributions were assumed for the 
velocity, temperature, and bubble density to ob- 
tain solutions for the steady-state case with 
given geometry and inlet conditions. Preliminary 
comparisons with experimental data obtained 
at the University of Eindhoven have been en- 
couraging. 

A number of existing two-phase dynamic 
models and data were examined from the view- 
point of the commonly adopted assumptions that 
(1) thermal equilibrium prevails and (2) local 
pressure variations can be neglected.”* Evi- 
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dence shows that failure to include nonequilib- 
rium and compressibility effects in the analysis 
can lead to erroneous results. Thus these as- 
sumptions clearly should be verified for future 
analyses of two-phase flow; the effect of com- 
pressibility and nonequilibrium, particularly, 
should be included in future analyses describ- 
ing the coolant behavior in nuciear reactors 
during rapid power excursions. 


PRESSURE-DROP PREDICTION 


In a survey lecture Levy’? discussed three 
types of analytical models for predicting pres- 
sure drop in two-phase flow: (1) one-dimen- 
Sional models, including the homogeneous and 
Martinelli—Nelson types and those based on 
momentum and energy exchange; (2) two-dimen- 
Sional models, illustrated by the Bankoff 
variable-density model and Levy’s mixing- 
length model; and (3) flow-pattern models, 
which are based on specified flow regimes 
(bubble, slug, annular, etc.). Levy concluded that 
considerably more research is needed on two- 
phase flow, especially in interface mechanics, 
and he singled out the flow-pattern model as the 
most promising approach for improving the pre- 
diction of pressure drop in two-phase flow. 


Boiling Heat Transfer 


The high heat fluxes obtained by boiling con- 
tinue to attract the attention of numerous in- 
vestigators. 


NUCLEATE BOILING 


Experimental and theoretical results were 
reported on the kinematic behavior of a bubble 
column originating at a single site during poil- 
ing of saturated liquid." From the study of 
vapor-removal kinematics during nucleate and 
film boiling, it was concluded that, during vapor 
removal by gravity, the volume flow rate is 
proportional to [g2(p, —9,)o]*8 in both boiling 
regimes, so long as the individual bubble column 
is undisturbed. 

The relation between local wall-temperature 
fluctuations and the formation of isolated bub- 
bles was discussed for pool-boiling of water at 
atmospheric pressure.”® Experiments indicate 
that the temperature drops correspond to the 
bubble growth and that the bubble departure has 
no perceptible effect on surface temperature. 
A proposed model for the bubble-formation 
mechanism was outlined. 
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From experiments and motion pictures con- 
cerning microconvection in boiling water, and 
from the determination of liquid volume ex- 
changed in bubble growth and vapor removal, 
it appears that in subcooled boiling the liquid 
acceleration is related to the mechanism of 
bubble condensation (not growth).”® In contrast 
to the situation with subcooled boiling, in satu- 
rated boiling there is no violent acceleration of 
liquid. Therefore only the buoyancy forces and 
the inertial forces related to the bubble growth 
period cause a Slow liquid transport. 


SIMULATED NUCLEATE AND FILM BOILING 


An experimental investigation simulated nu- 
cleate and film boiling by bubbling gas from a 
heated porous disk in a pool of water at atmo- 
spheric pressure.”’ From studies of how the 
heat-transfer coefficients vary with air-injec- 
tion rate, water temperature, diameter of the 
porous disk, and ratio of the pool diameter to 
the porous-disk diameter, it was concluded 
generally that: (1) an increase in the air-injec- 
tion rate is accompanied initially by anincrease 
in the heat-transfer coefficient; (2) when air 
injection reaches approximately 1.3 ft/sec, a 
transition from bubbles to a layer of air cccurs, 
and the heat-transfer coefficient decreases 
(see Fig. 7), simulating the transition from 
nucleate to film boiling; (3) with a constant air- 
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Fig. 7 Transition from nucleate to film boiling, sim- 
ulated by bubbling various amounts of air through a 
porous disk in a water pool, causes heat-transfer 
coefficient to decrease when air changes from bub- 
bles to a layer of air.?” For the test shown, disk tem- 
perature was 112°F and diameter was 2.5 in., pool 
diameter was 4.75 in., Prandtl number was 4.0, and 
Crassman number was 1.8 x 10", 
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injection rate, a rise in water temperature 
causes the heat-transfer coefficient to increase; 
and (4) variations in disk diameter or pool 
diameter do not significantly affect the heat- 
transfer coefficient. 

In a lecture on heat-transfer regimes, Ku- 
tateladze?® reviewed various aspects of his own 
criteria for predicting the transition from nu- 
cleate to film boiling. He made an interesting 
comparison between the actual boiling process 
and the one simulated by bubbling an inert gas 
through a heated surface that is transferring 
heat to a nonboiling liquid. In agreement with 
Akturk,”’ he concluded that the analogy between 
the two is good. Kutateladze discussed the 
phenomena associated with “burnout” or heat- 
transfer “crisis” under conditions of high- 
velocity liquid flow. He also presented useful 
dimensionless groups for defining two-phase 
flow patterns. 


FILM BOILING 


A study of interface dynamics in stable film 
boiling emphasized the frequency of generation 
of waves at a stagnation point.”® Analysis indi- 
cated that the phenomenon is dominated by ac- 
celeration effects in the vicinity of the stag- 
nation point, rather than by viscous effects. 
Photographic and motion-picture studies showed 
that the effect of the stagnation point is felt 
throughout the laminar-flow region. 

Another problem related to film boiling — 
that of waterdrops resting or impinging on a 
hot metal surface (Leidenfrost phenomenon) — 
was examined for measurement of: (1) evapora- 
tion rate of drops lying onahot surface, (2) heat 
transfer to droplets during their impact with a 
hot surface, and (3) heat transfer from a hot 
surface to a transverse spray flow.*” When the 
vapor layer that forms between the liquid and 
the hot wall is thicker than the roughness of the 
metal surface, the measured evaporation rate 
appears to be predictable by a simple model. 


Entrainment 


The problem of entrainment of a liquid ina 
gas, or vice versa, caused by incomplete separa- 
tion has interested the nuclear and chemical 
industries for some time. 

Because of concern with the detrimental effect 
of steam entrainment in the downcomer of a 
natural-circulation boiling reactor, data were 
taken in both a large-scale atmospheric-pres- 





oS ee a | 


\> p— =F 


a ee) ee oe eo 


=~ FF 


Winter 1965-1966 


sure air—water loop and in a high-pressure 
steam—water loop.*! The entrainment, repre- 
sented by the void volume in the downcomer, was 
given in terms of the independent variables: 
superficial downcomer water velocity, super- 
ficial riser steam velocity, water level, riser— 
downcomer area ratio, and pressure. Statis- 
tical correlation equations useful for reactor 
design were presented. 

The opposite problem of liquid entrainment 
in the vapor phase, which may occur when vapor 
is forced through a liquid region above an orifice 
of a sieve tray, was treated.* The mechanism 
by which entrainment is formed was reported, 
and a model based on this mechanism was form- 
ulated for predicting and extrapolating entrain- 
ment data. 

The maximum water flow possible, without 
gas entrainment in vertical pipes froma vessel, 
was graphed as a function of downcomer diam- 
eter.'° The experiments involved various gases 
and liquids and pipe diameters ranging from 1 
to 12 in. Figure 8 shows the effect of the pipe 
diameter on the critical flow rate. 
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Fig. 8 Effect of pipe diameter on critical flow rate 
for air—water mixture? 


Flooding 


When a liquid flows downward and gas flows 
upward in a vertical passage, there is a limit- 
ing condition known as “flooding.” Two papers 
dealt specifically with this problem. 

Experiments were conducted to study various 
aspects of the transition from falling-liquid-film 
flow to climbing-film flow in the presence of a 
rising gas stream.” The flcoding transition was 
found to be a strong function of the channel 
length. Figure 9 shows that the flooding veloci- 
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ties obtained were usually higher than those 
predicted by the theories of Schutt* or Davidson 
and Shearer® and that they spanned the pre- 
dicted values of Wallis.*® Continuous measure- 
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mospheric pressure depends on the test-section 
length and usually exceeds theoretical predictions.® 


ments of the liquid-film thickness and the 
pressure gradient resulted in adetailed charac- 
terization of the flooding phenomenon. 

Because of concern that flooding could be 
detrimental to the cooling of fuel elements in a 
nuclear reactor during emergencies, a series 
of experiments was conducted on single tubes 
and rod clusters.*’ Flooding was found to depend 
strongly on geometry. Although the nondimen- 
sional velocities proposed by Wallis® were use- 
ful in correlating the data, it was concluded that 
the mechanism of flooding is far from being well 
understood. 


Critical Flow 


Critical and high-speed two-phase flows are 
important in applications suchas nuclear reactor 
safety, in the analysis of flow from geothermal 
wells, and in devices for metering two-phase 
mixtures. 
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CHANNEL FLOW 


Experiments on critical air—water flows in 
tubes and channels at low pressures were re- 
ported for a large range of qualities.*® Data 
included critical flow rates, critical pressures, 
and exit void fractions. As Fig. 10 shows, the 
void-fraction measurements lead to much lower 
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Fig. 10 Experimental and calculated velocity ratios 
of air—water mixtures for 17-psia critical pres- 
Sure. 


values of the velocity ratios than had been pre- 
dicted by empirical and analytical expressions 
used in previous critical-flow models. The pre- 
vious models** *! agreed well with the previous 
data because the high values of velocity ratios 
compensate for metastability effects in the 
models used. A nonequilibrium-slip model was 
presented that agrees well with existing steam — 
water data. 

The vapor choking model was recommended 
as being applicable to thermal-equilibrium criti- 
cal mist flow, as in geothermal wells.*” Equa- 
tions were written suitable for a step-by-step 
solution of the mixture properties as the well 
exit is approached. Although exit slip and effec- 
tive flow areas were calculated for critical- 
flow conditions, no general test of the theory was 
possible because of the lack of suitable experi- 
mental data. 


NOZZLE FLOW 


The difference in mass-limiting two-phase 
flow in straight tubes and nozzles was discussed, 
with particular emphasis on the effect of back- 
pressure on the mass flow rate and pressure 
profiles when mass-limiting flow conditions are 
approached.** Figure 11 shows some of the data 
and indicates that each geometry has different 
characteristics at the near mass-limiting con- 
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Fig. 11 Mass flow as a function of exit-plane-to- 


veceiver pressure differential.43 Flow-control tap 
conditions were 53.8 lb;/sq in. and 0.116 quality; the 
vertical arrows show the approximate point of change 
of slope. 


ditions. Although the straight tube shows con- 
ventional choking characteristics, large pres- 
sure differentials between the exit plane and 
the receiver are required to produce mass- 
limiting conditions in nozzles. 


ORIFICE FLOW 


Data were presented for the flow of steam— 
water mixtures through sharp-edged orifices; 
correlations of the data are useful for flowmeter 
applications.“ For steam—water mixtures with 
steam qualities up to 10% and with 24 to 198 
psia upstream pressures of the orifice plates, 
the data can be correlated within 10% by: 


M,5/M, — 1=a(YF)?® 
where Y = (X/1 —X)(M,y/Myy) 


alg OUS n+i/n—-1 
F=1.414 ee (" . ") 
n 2 


a=1.6d°"4 
b=0.91 —log a 


Orifice inclination had no effect on the data. 


Two-Phase Flow Dynamics 


The transient behavior of two-phase flow sys- 
tems, which was considered in a number of 
papers, is of great importance in the design, 
control, and safety of boilers and nuclear power 
reactors, as well as in the performance of gas— 
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Fig. 12 Power-to-void transfer functions measured in steam—water mixtures within a vertical 
rectangular heated test section. For the three runs shown (left to right), average void fractions 
were 46.4, 34.5, and 52.6%, and the axial distances were 0.834, 0.467, and 0.834, respectively. 


liquid contacting operations such as distillation 
and gas absorption. 


NUCLEAR REACTOR APPLICATIONS 


Power-to-void transfer functions were mea- 
sured in steam—water mixtures within a ver- 
tical rectangular test section that was heated 
electrically;“* Fig. 12 shows typical data.” 

Transient-void data were obtained by a “one- 
shot” gamma-attenuation technique; conditions 
investigated included 200 to 800 psia pressures, 
2.52 to 3.78 ft/sec inlet flows, 2.28 to 9.12 x 
10‘ Btu/(hr)/(sq ft) heat fluxes, and0.8to12.4°F 
inlet subcoolings. Because the data deviated con- 
siderably from previous theoretical models, ‘’~*° 
a new model was constructed. 

Similar transfer-function measurements, but 
limited to 140 psia and to tubular test sections, 
were also reported.*! Transient-steam- volume 
fractions were obtained by a dilatation technique 
in closed-loop operation. Good agreement was 
claimed between the data and available models 
for two-phase flow dynamics. 

Consideration was given toneglecting body and 
friction forces in the mathematical description 
of unsteady steam —water flow in a heated chan- 
nel. Although the frequency responses of the 
void fraction and steam and water velocities 


were computed, a numerical example, based on 
the Halden Boiling Water Reactor, illustrated the 
need for including the body and friction forces in 
two-phase flow analysis. 

The onset of flow instability in a natural- 
circulation loop with two electrically heated 
parallel round ducts was discussed.*® Uniform 
power distribution was applied along the 4300- 
mm heated length of 20-mm-diameter channels 
at 50-bars pressure with 6°C subcooling. The 
data were discussed in terms of investigated 
loss coefficients in different parts of the loop. 

The dynamics of two-phase systems were 
examined in a study that concluded that flow 
oscillations cannot occur unless a maximum in 
the density—enthalpy product is obtained in the 
heater 4,55 


DISTILLATION APPLICATIONS 


In contrast to the papers that dealt primarily 
with nuclear reactor applications, several pa- 
pers considered the dynamic behavior of gas— 
liquid contacting devices as used in distillation 
operations. One paper described a two-dimen- 
sional hydrodynamic-stability analysis of the 
oscillations in a gas—liquid mixture on a sieve 
plate.*® The solution, which is based on numer- 
ous assumptions, indicated that neutral andam- 
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plified oscillations can occur only in distinct 
regions of the wave number. 

The need for the sieve-tray design to be based 
on dynamic, rather than on steady-state, con- 
cepts was illustrated by experiments that showed 
large fluctuations in froth pressure drop in 
normal sieve-tray operation.5” The amplitude 
of the fluctuations, which increased with the gas 
rate, reached 80% of the mean pressure drop. 

The liquid on the sieve tray was demonstrated 
to behave in either a stable or an unstable way, 
depending on the sieve-tray geometry and the 
resulting flow regimes. In the unstable regime 
the liquid-level wavelengths and oscillation fre- 
quencies were studied as functions of liquid 
height and gas flow. 

In view of the enormous complexity of the 
many two-phase flow phenomena, it seems ap- 
propriate to conclude that many of the sym- 
posium papers are preliminary and intermediate 
and that many more investigations are necessary 
to obtain data suitable for various design prob- 
lems. 


Symbols 


a, = sonic liquid velocity 
ay = sonic mixture velocity 
> = Crassman number 
= diameter of tube 
d = diameter of orifice 
G = mass flow 
g = gravitational constant 
k = velocity ratio =u¢/u, 
Mg = gas mass rate 
M, = liquid mass rate during two-phase flow 
M,n= liquid mass rate by weight during single-phase 
flow with unit contraction coefficient and two- 
phase pressure drop 
Mzs = liquid mass rate by weight during single-phase 
flow with normal contraction coefficient and 
two-phase pressure drop 
Myy= vapor mass rate by weight during single-phase 
flow with unit contraction coefficient and two- 
phase pressure drop 
n = index of isentropic expansion 
P = pressure 
Po = critical pressure 
y = ratio of backpressure to initial pressure 
Reg = gas Reynolds number 
Re, = liquid Reynolds number 
S = dimensionless distance from wall 
= average gas velocity 
T; = porous surface temperature 
Vsq = superficial gas velocity 
, = superficial liquid velocity 
Np, = Prandtl number 
We = Weber number 
X = quality 
q@ = void fraction 
£ = volumetric gas-flow fraction 
= gas density 
p,, = liquid density 


> 
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& = dimensionless axial distance 
Ug = gas viscosity 
uy = liquid viscosity 

o = surface tension 
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How Coolant Dynamics Affect Fast Reactor Safety 
By Ralph M. Singer 


The safe operation of a fast reactor can be 
affected adversely by two general types of 
occurrences: (1) loss of cooling capability and 
(2) power excursions caused by reactivity in- 
sertions. The good safety experience with fast 
reactors is based on substantial design and 
operations-management efforts to assure re- 
tention of cooling capability and to prevent 
unplanned insertions of reactivity. However, 
as more and larger fast reactors are planned, 
it becomes increasingly desirable to under- 
stand and predict the consequences of such 
events if they should occur despite the pre- 
cautions. Then inherent design features can be 
incorporated to assure reactor safety. 

If the extent of the consequences is to be 
understood and predicted when cooling capacity 
is lost or unwanted reactivity is inserted, it is 
necessary to understand and predict the role of 
each mechanism that would affect reactor be- 
havior —and also to know the precise sequence 
of these various mechanisms as well as the 
effects of their interactions. Then it would be 
possible to predict which specific accidents, 
such as pump failures, fuel-rod bowing, or 
coolant-system leaks, could damage the system. 

Although the probable sequence of events that 
would follow cooling loss or a power excursion 
in a fast reactor is of vital interest to reactor 
designers, little information is available on 
specific mechanisms relating to coolant dy- 
namics and heat transfer. 


Sequence of Events 


Figure 1 shows the possible sequences of 
events that would follow a loss of cooling or a 
power excursion in a fast reactor. Although 
many Sequences are possible, the sequence that 
actually would occur would depend on the 
specific details of the accident and on the re- 
actor system itself. 


EXAMPLE 


A sufficiently rapid rate of insertion of re- 
activity can cause the fuel and cladding tem- 
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Fig. 1 Possible sequence of events that could follow 
loss of cooling or a reactivity-insertion-induced 
power surge ina fast reactor, excluding veactivity- 
feedback effects. 


peratures to rise so fast and high that the fuel 
and/or cladding will melt before enough energy 
transfers to the coolant to cause boiling. The 
melting of the fuel can result in at least three 
events: (1) the coolant may vaporize, (2) a 
pressure pulse may occur, or (3) the core may 
disassemble. If the coolant vaporizes after the 
core melts, there may result a pressure pulse 
and coolant expulsion, either of which can re- 
sult in additional core damage. Expulsion of 
coolant could cause an additional reactivity 
increase that could lead to further damage. 

If the initial reactivity insertion rate is 
moderate, sufficient energy could be trans- 
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ferred to the coolant to vaporize it before the 
fuel melts; then the vapor could insulate the 
fuel, causing it to melt. Alternatively, a tran- 
sient power surge in the core could suddenly 
expand the fuel and pressurize the coolant so 
that the coolant is directly expelled, which 
could cause the core to melt and disassemble. 

Other sequences of events are undoubtedly 
possible, but the example indicates the com- 
plexity of the situation and the need for under- 
standing not only the individual processes but 
also the interrelations among them. Itis worth- 
while to consider the status of our knowledge 
of transient heating of the core and single- 
phase behavior of the coolant, transient heat 
transfer and coolant hydrodynamics, transient 
vaporization of the coolant, pressure-pulse 
generation, and coolant expulsion. 


Transient Core Heating and Single- 
Phase Behavior 


If the energy generated in a reactor fuel is 
to be transferred to the reactor coolant, the 
heat must be conducted through the fuel ma- 
terial, the bonding material (if any), and the 
cladding material. Because the space-time 
variation of the temperature in the core must 
be known if core behavior is to be predicted, 
the first step in a safety analysis should be a 
study of transient heat conduction in multi- 
dimensional nonhomogeneous media, including 
the effects of internal energy generation and 
phase changes. 

An excellent survey has been made of the 
available digital-computer codes which can 
calculate transient temperatures in various 
geometric regions and which consider coupling 
to reactor kinetics and coolant hydrodynamics.’ 
The survey emphasizes the four codes, ARGUS, 
FORE, TIGER-V, and TER-4, and discusses 
their applications and limitations in detail. An 
additional code, TRANSFUGUE-I, deals with 
transient hydrodynamics and heat transfer.’ 
All existing codes involve approximations and 
limitations that restrict their general applica- 
bility. Perhaps the most glaring deficiency in 
all the codes is their complete neglect of 
compressibility effects (critical flow and low 
sonic velocity) in two-phase mixtures. In cer- 
tain situations sonic velocity in two-phase 
mixtures can be as low as 5 ft/sec, which can 
cause channel choking as soon as voids are 
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formed. A recent paper discusses this prob- 
lem analytically.° 

Another recent paper, which presents a new 
numerical method for solving multiregion con- 
duction problems, claims to combine the com- 
putational ease of explicit methods with the 
stability of implicit methods.’ Use of this tech- 
nique should be investigated in any future 
transient analysis codes. A reasonably detailed 
survey’ of methods for predicting fuel-element 
and coolant behavior under reactor transient 
conditions summarizes studies conducted be- 
fore 1961. 

The prediction of incipient fuel melting or 
coolant vaporization has been simplified by the 
development of methods for predicting tem- 
perature transients of the fuel and single-phase 
coolant.® The simplified analytical techniques 
can be used to describe heat transfer in typical 
fuel channels and to obtain a rational extrap- 
olation of existing empirical information to 
boiling sodium flow through heated channels. 


Coolant Vaporization 


When the space-time history of tempera- 
tures in a reactor core has been predicted for 
a specific incident, the next step is to deter- 
mine whether any phase changes occur in the 
fuel or coolant. Prediction of the conditions 
under which the coolant will vaporize, and thus 
prediction of void-reactivity and other effects, 
requires a knowledge of (1) the amount of 
superheat needed to initiate boiling and (2) the 
factors that affect the superheat. For liquid 
metals, large superheats may be necessary to 
initiate boiling —and, when voids start to form, 
void growth can be extremely rapid." Al- 
though work that was conducted before 1965 has 
been reviewed,’ only a very limited amount of 
experimental information is available. 

Once nucleate boiling begins, liquid metals 
have been observed to behave very unstably, 
with large variations in wall temperature and 
fluctuation between nonboiling and boiling 
states.*.!0 Apparently this instability is caused 
by the physical characteristics of the heated 
surface as well as by the properties of the 
fluid.!° However, this instability has not been 
observed by all experimenters. 

Data on transient boiling (pool and forced 
flow) of liquid metals are vital to designers 
for prediction of the space-time history of 
voids and temperatures. Although there are no 
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data for liquid-metal systems undergoing rapid 
heating transients, it might be possible to 
extrapolate to liquid metals the limited amount 
of information that exists on the transient boil- 
ing of water from a ribbon.!!~!® Extrapolation, 
of course, has inherent uncertainties and should 
be pursued cautiously. 

A recent analytical description of transient 
boiling during fast. reactor power transients 
leads to the conclusion that: (1) voiding in a 
reactor channel starts near the upper boundary 
of the core and grows axially along the fuel 
channel, (2) vapor production is so rapid that 
transient self-pressurization of the fuel channel 
occurs, and (3) the resulting two-phase flow is 
strongly unstable hydrodynamically.'? Under 
certain circumstances, power sufficient just to 
initiate boiling can result in complete voiding 
of the channel without further increase in 
power. 

The importance of coolant vaporization in 
fuel-meltdown studies, suggestions for needed 
work, and some answers for the following 
questions are described”” in a survey report: 
(1) What abnormal conditions may lead to 
coolant vaporization? (2) Which phenomena as- 
sociated with coolant vaporization are most 
important in each of the various possible ac- 
cidents ? (3) What are the important character- 
istics of each phenomenon, and what methods 
are available to predict them ? 


Generation of Pressure Pulses 


For a safe design to be obtained, it is es- 
sential to know what pressures are attainable 
if the reactor operates under abnormal con- 
ditions. The SPERT tests show that destructive 
pressure transients can be generated in water 
reactors, although the mechanisms Causing 
them are still matters for discussion. Pro- 
grams that are under way at Argonne National 
Laboratory,” TRW Space Technology Labora- 
tories,” and elsewhere are aimed at deter- 
mining and understanding these: mechanisms. 
Although experimental work so far has been 
limited to water systems, many of the pressure- 
generation mechanisms considered could apply 
equally well to liquid metals. Several of these 
mechanisms involve: sudden thermal expansion 
of the fuel and immediately adjacent coolant, 
rapid water —metal reactions that evolve hydro- 
gen, pressure amplification caused by the col- 
lapsing of voids, and sudden void formation. 
However, much experimental and theoretical 
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work must still be done before it is possible to 
predict accurately the severity of pressure 
transients that can result from specific ac- 
cidents in specific reactors.. 


Coolant Expulsion 


After a pressure pulse is generated in the 
coolant, the coolant will be expelled from the 
vicinity of the high-pressure source. Before 
the expulsion velocity can be predicted (as- 
suming that the pressure-generating mechanism 
is known), it is necessary to understand the 
transient mechanics of two-phase flow. Al- 
though the theoretical tools are not yet devel- 
oped sufficiently for two-phase flow, the ex- 
pulsion velocity of a single-phase fluid can be 
estimated if enough is known about the initial 
pressure pulse. It is not known whether such 
an estimate would be accurate enough for 
confident predictions of reactivity feedback 
effects. 

Apparently there is no experimental infor- 
mation on liquid-metal expulsion and only ex- 
tremely limited data on water expulsion.”»*4 
Motion pictures conclusively show that critical 
flow (choking) occurred in a two-phase system 
in which superheated water was ejected from a 
vertical tube."4 The water was superheated by 
means of overpressurizing; the pressure sud- 
denly decreased by rupturing a diaphragm, and 
the result was boiling and ejection of the water. 
An expulsion experiment under way at Argonne 
employs a large power surge in an electric 
heater to boil water in an annulus. It will 
produce data on expulsion velocities, pressure 
generation, and void growth rates and will lead 
into related experiments with sodium. 


Conclusions 


From this brief review of the effects of 
coolant dynamics on fast-reactor safety, it can 
be concluded that: 


e The transient thermal and hydrodynamic 
behavior of the core and coolant is reasonably 
understood if no phase changes occur. 

e There is a lack of data on how fuel melt- 
ing (with or without cladding failure) affects the 
thermal and hydrodynamic behavior of the 
coolant. 

e Data are needed on the degree of super- 
heat needed to cause nucleate boiling of liquid 
metals under various conditions. 
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e Data on the transient nucleate boiling of 
water from a ribbon have not been rationally 
extrapolated to predict the transient nucleate 
boiling of liquid metals from a fuel surface (if 
at all possible). 

e Data are lacking for the transient boiling 
(pool and forced flow) of liquid metals after a 
power surge; data on heat-transfer rates and 
void-formation rates are required. 

eA generally accepted explanation of the 
mechanisms that can generate destructive pres- 
sure transients is not available; the available 
data for water systems should be rationally 
extrapolated to liquid-metal systems, and data 
should be obtained on liquid-metal systems. 

e Data on coolant expulsion rates are ex- 
tremely sparse. 


The transient behavior of the coolant is 
probably the least understood phenomenon oc- 
curring in a reactor. However, the papers and 
discussions on the subject at the Conference on 
Safety, Fuels, and Core Design in Large Fast 
Power Reactors, which was held at Argonne in 
October 1965, indicate that there was general 
agreement among U.S. and European engi- 
neers concerning the proper direction of re- 
search needed to resolve the problems. Thus 
the results of current and planned research 
programs soon will yield an improved under- 
standing of these complex problems. 
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UHTREX Operation Near 


By Barton M. Hoglund 


A novel gas-cooled reactor concept, with the 
acronym UHTREX (Ultra High Temperature 
Reactor Experiment), is expected to achieve 
criticality in October 1966 at the Los Alamos 
Scientific Laboratory. Among its unique fea- 
tures of interest to reactor designers are: 

e An annular core that is rotatable to facili- 
tate refueling during reactor operation, 

e Articulated control rods that can be in- 
serted directly into the core regardless of its 
angular position, 

@ Simple graphite fuel elements that hold 
promise for fuel economies, but that lead to 

e Normal operation with a contaminated pri- 
mary coolant, and 

® Use of an extensive on-line coolant- cleanup 
system to limit the amount of contamination in 
the coolant. 


Purpose and Potential 


The UHTREX, a 3-Mw(t) helium-cooled graph- 
ite-moderated reactor experiment, has been 
built so that operating experience can be ac- 
quired with a gas-cooled reactor system that 
employs simple graphite type fuel elements of 
various types. In planning the experimental pro- 
gram to be conducted in UHTREX, Los Alamos 
is considering using the fuel-loading flexibility 
of the reactor to determine the operational per- 
formance of various fuel materials. For ex- 
ample, to determine how neutron economy would 
benefit if gaseous fission products were re- 
leased from the fuel during reactor operation, 
the experimental program might employ fuel 
elements that release fission products more 
rapidly than coated-particle fuels that will be 
uSed initially. 
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“Should fuel simplification outweigh the diffi- 
culties arising from the operation of a con- 
taminated coolant circuit, UHTREX may provide 
the basis for significant economic gains to be 
made in future gas-cooled reactors for both 
power and process-heat applications.” In addi- 
tion to power and high-temperature process 
heat for such uses as coal gasification, this 
concept might be a useful energy source for 
magnetohydrodynamic (MHD) devices or for 
high-temperature Brayton-cycle applications. 
The hope has been expressed! that, by develop- 
ment of high-temperature reactor technology 
beyond current applications, the applications 
technology would be stimulated to “catch up.” 


HISTORY 


UHTREX is a continuation of the TURRET re- 
actor experiment,’ which was first proposed as 
a nitrogen-cooled reactor to drive a gas turbine 
in a closed cycle.® The proposal later was re- 
vised‘ to make TURRET a helium-cooled reac- 
tor having an operating pressure of 500 psia and 
an exit temperature of 2400°F. The TURRET 
program was terminated in January 1961; the 
UHTREX program began in May 1961. The 
UHTREX has been constructed on a site about 
2.5 miles southeast of Los Alamos. 


Design Description 


The UHTREX facility (Fig. 1) includes oper- 
ating and control areas, a fuel-handling area, 
and a containment structure. All pressurized 
equipment containing contaminated helium is 
within the containment, which is a single con- 
crete structure that is lined with a steel mem- 
brane and is divided into (1) a separately 
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Fig. 1 


The UHTREX facility consists of reactor primary system, containment structure, and con- 


trol and service equipment rooms. The containment encloses reactor cell, primary coolant loop, 
fuel-loading apparatus, and fuel-discharge and gas-cleanup rooms.° 


Shielded reactor cell for the reactor, primary 
coolant loop, and fuel-loading apparatus; (2) a 
fuel discharge room; (3) a gas-cleanup room; 
and (4) two gas-cleanup piping rooms. 

The UHTREX reactor (Figs. 2 and 3), which 
has a vertical cylindrical core surrounded by a 
reflector and insulation, is housed in a spheri- 
cal pressure vessel with top-mounted control- 
rod drives. Fuel elements are loaded through 
one side of the reactor and discharged from the 
bottom. Table 1 (page 5) summarizes the 
parameters. 


Fuel Elements 


Each of the initial fuel elements to be used in 
UHTREX consists of an extruded hollow cylinder 
made of a mixture of graphite and93% enriched 
UC, particles. The 5.5-in.-long cylinder has an 
inside diameter of 0.5 in. and an outside diam- 
eter of 1 in. Four of these fuel cylinders can be 
positioned end-to-end in the reactor horizontal 
fuel channels, each of which is 23.5 in. long and 
has a diameter of 1.1 in. 


COATED-PARTICLE FUEL 

The UC, particles are coated with pyrolytic 
carbon (PyC) to retard the release of fission 
products. The diameter range of the coated 


particles is about 350 to 450 yu, resulting from 
100 to 110 pof PyC coating on 147- to 208-yu 
UC, cores. This range represents a compromise 
between (1) the requirements for good particle 
strength and fission-product retention and (2) 
the need to keep the particle size and volume 
fraction of particles small to avoid difficulties 
in fabrication by the extrusion process. The 
PyC on the UC, particles is a General Atomic 
“Triplex” coating in which the 20- to 30-yinner 
coating is a buffer layer of low-density PyC to 
catch fission recoils, to dissipate stresses 
caused by differential dimensional changes, and 
to provide a void volume that will minimize 
fission-gas pressure buildup. The 35- to 45-y 
middle layer consists of “isotropic” (randomly 
oriented) PyC with low permeability to fission 
gases and good resistance to radiation damage. 
The outer layer is fine-grained columnar PyC. 
The interface between the isotropic and co- 
lumnar layers is intended to retard the propaga- 
tion of any cracks that might develop.® Fuel- 
element properties and fabrication methods and 
equipment have been described.’ 


MEASURING URANIUM CONTENT 


The uranium content of each fuel element is 
determined nondestructively after fabrication 
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Fuel discharge 


Fig. 2 Vertical cutaway of reactor shows horizontal fuel channels in annular core and vertical slot 
in stationary tapered central plug, through which discharged fuel elements drop when displayed by 


new element inserted by vams at right.'9 
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by measurement of the gamma emission from der at 15° intervals in 13 horizontal planes, are 

the uranium with a NaI scintillator detector that separated vertically 3 in. center-to-center. 

completely surrounds the element. The inte- Surrounding the core is a 16-in.-thick reflector 
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Fig. 3 Horizontal cutaway shows stationary central plug with fuel-discharge slot and eight channels 
for regulating rods (plus center channel for startup source), rotatable annular core with radial fuel 
channels and 24 vertical shutdown-rod channels, the outer reflector, insulation, and vessel wall.'9 


grated current produced in scintillator photo- made of 4 in. of graphite and 12 in. of dense un- 
tubes relates directly to the uranium content. graphitized carbon. Porous carbon between the 
Less than 2 min is required for a measurement reflector and the 13.5-ft-OD carbon-steel pres- 
with a relative standard deviation® of 0.3%. sure vessel serves as insulation. 


The hollow axial region of the core containsa 


Reactor Core tapered graphite plug; the taper leaves a plenum 


The annular rotatable core is a hollow verti- between the plug and the core through which the 
cal cylinder (70 in. in outside diameter by 23 in. helium coolant enters the core to flow radially 
in inside diameter by 39 in. high) in which outward through the fuel channels. The entire 
graphite serves as both the structural material annular core, together with portions of the re- 
and the moderator. The cylinder is made up of flector and insulation, can be rotated about the 
five concentric hollow cylinders so as to limit cylindrical axis so that fuel can be loaded and 
the thermal stresses across each piece of unloaded from the various channels. The mov- 
graphite. The fuel channels, which are 1.1-in.- able portions (Fig. 4) rotate ona large-diameter 


diameter holes bored radially through the cylin- dry-film-lubricated ball bearing whose integral 
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Table 1 SUMMARY OF UHTREX REACTOR CHARACTERISTICS!:? 





Reactor type 
Rated power 
Fuel elements 


Moderator 
Reflector 
Reactor vessel 
Fuel channels 


Critical loading 


Temperature coefficient 
of reactivity 


Power coefficient 
of reactivity 

Thermal-neutron flux 

Fuel specific power 

Core power density 

Fuel burnup 

Refueling 

Rate 

Excess reactivity required 
to compensate for tem- 
perature coefficient 

Control-rods material 


Total worth 
Withdrawal 

Coolant 
maximum temperature 


Experimental, high-temperature, gas-cooled, unclad fuel 

3 Mwit) 

Graphite cylinders, 1 in. OD, 0.5 in. ID, 5.5 in. long; pyrocarbon- 
coated fully enriched UC, particles dispersed in graphite 

Vertical graphite cylinder, 39 in. high, 70 in. OD, 23 in. ID 

4 in. of graphite plus 12 in. of dense carbon 

Carbon-steel sphere, 13.2 ft ID, 1.75-in. wall, internally insulated 

312 radial, horizontal holes, each of 1.1 in. diameter, 23.5 in. 
length; 4 fuel elements/channel 

11 kg 93% enriched uranium at design power and temperature and 
equilibrium poison; 5.7 kg 93% enriched uranium at zero power 
and 70°F 

—1.52 x 107 Ak/°C in operating region (prompt) 

—-1.05 x 10 Ak/°C in operating region 

—1.2 x 10 Ak/°C average from ambient to operating temperature 

—9,2 x 10 Ak/k/kw (-0.13 ¢/kw) 


2 to 4 x 10'3 neutrons/(cm*)(sec) 

272 watts/g of 755U (average) 

1.3 watts/cm 

10 to 50% 

At power 

1 to 6 fuel elements/day; cycle time about 6 min/element 
~20% Ak 


B,C in central rods; 1.9% boron- stainless steel in articulated core 
rods 

22% Ak 

18 in./min 

Helium at 500 psia max.; 10,250 lb/hr max. 

1600°F inlet; 2400°F outlet 
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ring gear is mounted at the bottom of the pres- 
sure vessel. 


Reactor Vessel 


The spherical pressure vessel supports the 
core, moderator, reflector, insulation, and ex- 
ternal components, including the fuel-loading 
mechanism and the control-rod drives (Fig. 2). 
The ASTM A-212B carbon-steel vessel is 13.4 
ft in outside diameter with a 1.75-in. wall; it is 
designed for operating conditions of 650°F and 
550 psig and for a test pressure of 825 psig. 

As shown in Figs. 2 and 3, penetrations for 
the fuel loader and helium outlet nozzle are on 
the center line at opposite sides of the vessel. 
The fuel-elevator housing and a massive sec- 
tion penetrated by fuel-loading holes, both inte- 
gral parts of the vessel wall, set up severe 
thermal stresses in the shell that are accommo- 
dated by a 3-in.-thick segment of the vessel 
wall. On the opposite side of the vessel is a 
short, 38-in.-ID nozzle for joining the recuper- 
ator to the pressure vessel. The vessel top is 
closed by an 11.7-ft-OD flanged cap that is 
bolted and seal welded; no gasket is used in the 
joint. The control-rod nozzles penetrate the 
cap. At the vessel bottom is a cylindrical ex- 


tension that has a side nozzle to the helium in- 
let line and a dished bottom head that is pene- 
trated by the fuel-discharge channel. There are 
two additional penetrations in the lower half of 
the sphere—one for instrumentation leads and 
one for the core indexing mechanism. 


Reactor Control 


Critical experiments show that the control 
system has a total rod worth of about 22% Ak. 
This is achieved by the insertion of 4 of the 12 
control rods into the rotatable core region. De- 
velopment of articulated control rods cleverly 
solved this alignment design problem. 

As shown in Fig. 3, the 12 control rods hang 
in vertical holes arranged in two circles cen- 
tered about the core axis;'’! the inner 8 holes 
are at a radius of 7.9 in. from the axis (and are 
in the stationary central plug), and the outer 24 
holes are 23 in. from the center (and thus are 
in the rotatable core). 


REGULATING RODS 


The inner eight control rods, whose channels 
are within the stationary central graphite plug, 
are calculated” to be worth 14% Ak. These cen- 
tral rods will be used for control and safety 
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Fig. 4 Annular core as well as upper and lower re- 
flectors and support block are mounted on support 
bearing that can be rotated to align horizontal fuel 
channels with refueling rams.?° 


during reactor operation as well as for shutting 
down the reactor. 

Because these central regulating rods are in 
the hot reactor atmosphere during operation, 
they are made of refractory materials (Fig. 5). 
Each regulating-rod assembly, whose diameter 
is 1 in. and active length is 40 in., consists of 
pressed B,C cylinders mounted on a 0,.25-in.- 
diameter niobium support rod that has a flange 
swaged onto its lower end (because niobium- 
weld areas embrittle when exposed to impure 
2400°F helium,® the support rod and end flange 
had to be swaged rather than welded). The 
regulating rod in Fig. 5 has a graphite sheath, 
but scram deceleration stresses have required 
that the design be changed to aniobium sheath.’® 
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Niobium support rod 
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7 Integral Nb flange 
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Fig. 5 The regulating rod (left), which moves within 
stationary central graphite region during reactor op- 
eration, is made of vefractory materials; in newer 
design, outer sheath is niobium rather than graphite. 
The shutdown rod (right, but not to same scale), which 
is made of type 304 stainless steel—1.9% boron alloy 
and is inserted into channels in rotatable core only 
during extended shutdowns or for emergencies, con- 
sists of 4-in.-long ball-jointed segments that can go 
into channel if rod and channel are aligned! within 


132". 
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CORE SHUTDOWN RODS 


The outer set of four control rods, whose 24 
channels are within the rotatable core region, 
normally will be completely withdrawn from the 
core during reactor operation—they will be in- 
serted only in case of an emergency or for ex- 
tended shutdowns, when the reactor cools enough 
that additional shutdown reactivity is required 
to override the negative temperature coefficient 
of reactivity. 

Analysis indicates that the eight central (in 
the stationary plug) regulating rods are ade- 
quate to bring the reactor subcritical at tem- 
peratures well below 900°F. Thus, by allowing 
the core to cool below 1800°F before inserting 
the core shutdown rods, the designers avoid 
having to make these articulated rods of re- 
fractory materials. 

Each of the core rods, made of type 304 
stainless steel—1.9% boron alloy, has a 1.4-in. 
outside diameter and a 40-in. active length.'* 
These four rods are made of assemblies of 4- 
in.-long ball-jointed segments (Fig. 5) to pro- 
vide flexibility so they can be inserted into the 
core even if the channels are not aligned with 
the rod-drive mechanisms. The core has 24 
vertical channels into which the four shutdown 
rods can be inserted (Fig. 3) to permit up to 
7.5° of rod-to-channel misalignment; these 
vertical channels are spaced at 15° intervals 
and pass between the horizontal radial fuel 
channels. If necessary, the core can be rotated 
as much as 15° with the articulated rods in- 
serted.'° Thermal-cycling tests show that all 
the ball joints remain free-acting throughout 
1542 cycles,!>»18 The maximum temperature at 
which these rods can be inserted into the core 
is 1800°F, above which significant socket scal- 
ing occurs. Contamination of the graphite mod- 
erator by boron from the shutdown rods will be 
negligible because core temperature will drop 
after shutdown and these articulated rods will 
be in only line contact with the holes. 


CONTROL-ROD DRIVES 

One basic type of control-rod-drive mecha- 
nism has been designed for both the central 
regulating rods and the articulated shutdown 
rods.’ The drive consists of an a-c hysteresis 
synchronous motor, with a d-c brake, driving a 
rack-and-pinion through a 940: 1 ratio planetary 
gear train at a single drive speed of 18 in./min. 
The gear rack is machined on a ball spline that 
is used for the vertical guide. The entire drive 
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mechanism is within the high-pressure primary 
coolant region, which eliminates seal problems. 

Free-fall scram of the eight central rods is 
triggered by interrupting electric power to a 
stationary magnetic clutch that normally en- 
gages the gear train to the pinion. The free- 
falling control rods are decelerated in the 
control-rod dashpot by '/.-in.-diameter hard- 
ened and ground steel balls that collectively 
act as a “semi-fluid.” Successful tests were 
run on a device with a 0.95-in.-diameter knob 
on the end of a 0.75-in.-diameter shatt that fell 
through a 1.625-in.-OD tube containing these 
balls.!" A “spring” formed by a 12-turn helical 
path cut into the wall of the impact tube eases 
the initial deceleration shock.'® 

Analytical results indicate that the reactor is 
stable and that the overall system is well be- 
haved. Typical thermal time constants are in 
the range® 100 to 1000 sec. Although thermal 
excursions in the core might arise from diffi- 
culties in the coolant loops, metallic compo- 
nents in the loops would be affected much more 
rapidly and radically than the refractory core 
would be. Thus action must be taken to protect 
the coolant loop from thermal effects well be- 
fore the reactor itself would suffer damage.’ 


Fuel-Handling System 


Figure 6 shows the entire fueling system for 
UHTREX. The short horizontal conveyor brings 
fuel elements from the fuel-handling cell; a 
dual-valve gas lock admits them into the con- 
tainment. The sloping conveyor raises each 
element to the top of the gas-lock sequence, 
through which it enters the reactor coolant 
system. The top ram in the loader mechanism 
advances the element into the fuel elevator on 
the side of the reactor vessel; the elevator 
positions it at any of the 13 levels of horizontal 
fuel channels. Then the electrically powered 
loading ram at the desired level pushes the ele- 
ment into the outer edge of the core and, in so 
doing, pushes the innermost element so that it 
drops through the vertical slot in the central 
core plug, passes down through the vessel ex- 
tension and lower gas lock, and is removed by 
the discharge conveyor. After the discharged 
element drops through the two valves to leave 
the containment enclosure, the long horizontal 
conveyor carries it into the fuel-handling cell. 
To align various fuel channels for refueling, the 
core can be rotated about its vertical axis by 
the electrically powered indexing mechanism."° 
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Fig. 6 Incoming fuel on conveyors enters containment and reactor system through two sets of gas- 
valve locks. The elevator (on right side of reactor vessel) aligns fuel element vertically with proper 
channel in core. When ram inserts new element, displaced spent element falls through central slot, 


leaves reactor,and is conveyed from primary system and containment through two sets of gas-valve 
locks.1° 
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GAS VALVES 


Considerable development work and testing 
have gone into the design of ball valves for the 
gas locks in the reactor coolant system. It was 
concluded’® that “it is not possible to devise an 
all-metallic seal with a leakage consistently 
less than 107° standard cm*/sec at 500 psi dif- 
ferential pressure under severe service condi- 
tions in which abrasive particles are introduced 
into the valve as it is operated through several 
thousand cycles. Soft metal seals were found to 
be effective initially, but the introduction of 
traces of particulate matter caused them to leak 
badly, even after cleaning wipers were em- 
ployed. The hard metal seals were not so ad- 
versely affected by abrasive particles but un- 
der 500 psi pressure differential, they never 
seal very well, even initially.” 

The prototype ball valve has a hard metal 
seal that holds pressure from either direction. 
Under a continuing series of endurance tests,'‘ 
this valve survived 2500 cycles under a 250-psi 
pressure differential with leakage of only 10~° 
em?/sec. After the seat was relapped, the valve 
endured 4500 more cycles with a 500-psi differ- 


PRIMARY LOOP 
10,250 tb/hr 
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ential before helium leakage exceeded 1.3 x 
10 cm*/sec. Primary coolant leakage past a 
primary valve is bled to a low-pressure re- 
ceiver for storage until it is pressurized and 
returned to the primary coolant system.! 


Cooling System 


At full reactor power, two closed-cycle he- 
lium loops in series (Fig. 7) remove 3 Mw of 
heat from the reactor core and dissipate it to 
the atmosphere. The primary coolant loop 
forms a closed system inside the containment. 
Heat passes from the primary loop to the sec- 
ondary cooling loop through the main heat ex- 
changer; the secondary loop takes the heat toan 
air-cooled heat exchanger outside the reactor 
building. Helium flow rates in the two loops are 
matched. 

In the primary coolant loop, 1600°F helium 
enters the reactor central plenum, from which 
it passes radially outward through the fuel 
channels, where it is heated to 2400°F. The he- 
lium is cooled to 1400°F as it passes through 
the hot side of the recuperator, one end of 
which is bolted and seal-welded directly to the 
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Fig. 7 In primary loop, hot helium from reactor passes through recuperator (which cools it 1000°F) 
and heat exchanger (which transfers heat to clean helium in secondary loop); blower circulates 
primary-system helium back to reactor through recuperator. The secondary clean-helium loop dis- 
sipates heat through radiator to atmosphere. Cleanup loop removes fission products from primary- 


system helium.’ 
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reactor vessel; the other end is connected to the 
main heat exchanger. 


THE RECUPERATOR 

This internally insulated regenerative heat 
exchanger made of graphite in a carbon-steel 
pressure shell cools the outlet helium to tem- 
peratures tolerable in the stainless-steel heat 
exchanger. Heat exchange in the recuperator is 
achieved in a 2- by 2- by 8.7-ft block of graphite 
that is penetrated by 440 holes bored parallel to 
the long axis. The 0.5-in.-diameter holes are 
drilled on '/,-in. centers arranged in 22 col- 
umns. The hot and cold gases flow countercur- 
rently through adjacent holes.’® 


MAIN HEAT EXCHANGER 


From the recuperator the helium flows to the 
main heat exchanger where it is cooled to 600°F 
by transferring heat to the secondary coolant 
loop. The main heat exchanger is made of type 
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316 stainless steel and has bayonet tubes to 
avoid thermal-stress problems. The 87 bayonet- 
tube assemblies each consist of three concentric 
tubes. When the 1400°F primary coolant enters 
the exchanger, it flows down the innermost tube, 
reverses direction at the bottom, and then flows 
up through the outer tube, where it exchanges 
heat with the secondary coolant in the main 
shell. The intermediate tube in the bayonet as- 
sembly serves only to provide astagnant helium 
gap that acts as a thermal insulator between the 
entering and exiting primary coolant; a sample 
wire that can be inserted into the annulus of the 
bayonet assembly can be withdrawn periodically 
for radiochemical analyses to measure the rate 
of fission-product buildup in the heat exchanger. 


MAIN BLOWER 

When the primary coolant leaves the heat ex- 
changer, it passes through the main circulating 
blower, then back to the recuperator inlet. The 
helium blower (Fig. 8) is driven by a three- 
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Fig. 8 The main blower raises pressure of 10,250 lb/hr of 510-psi helium 8 psi between heat ex- 
changer and recuperator inlet. Helium serves as lubricant for blower bearings, and a water-cooled 


heat exchanger (not shown) cools by-passed helium that in turn cools shaft journals and motor. 
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phase two-pole squirrel-cage induction motor fission-product contamination in the primary 
rated at 65 hp at 1200 rpm. The open-face im- loop is reduced, by the continuous purification 
peller on the centrifugal blower has an 11.6-in. of 102 lb of the helium stream per hour in the 
diameter and 19 vanes. Helium-lubricated hy- gas-cleanup system, a sidestream of the pri- 
drodynamic bearings support the rotor, and hy- mary loop (Fig. 9). 


drostatic pressure is applied during starts and 
stops; a centrifugal fan circulates helium through 
a water-cooled heat exchanger to cool the shaft 
journals and motor.’ 4 variable-frequency al- 
ternator, magnetically coupled to two 75-hp 60- 
cycle motors, drives the blower motor. Each 
75-hp motor can carry the alternator load and 
is supplied by a separate power source to pre- 
vent loss of pumping power if one of the main 
power feeders fails. With helium at 500 psi and 
600°F, the blower is designed to deliver 10,250 
lb/hr with a pressure rise of 8 psi. 


Gas-Cleanup System 


Helium from the primary loop enters the 
cleanup system through a replaceable stainless- 
steel filter that removes particulate matter as 
small as 0.5» in diameter, is heated to 1200°F 
in an electrical preheater, then flows into a 
copper oxide bed, where carbon and hydrogen 
compounds in the helium are oxidized to CO, 
and H,O, respectively. After leaving the copper 
oxide bed, the helium is cooled to 100°F and 
enters the tubes of a water-cooled decay heat 
exchanger. Activated carbon packed in these 
tubes delays passage of the krypton and xenon 
fission products by gas chromatography and ad- 
sorbs any iodine or bromine fission products. 





The chemical inertness of the coolant is pre- The radioactive decay of the fixed halogens and 
served, and the amount of volatile long-life delayed noble gases liberate about 12 kwof heat. 
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Fig. 9 Gas-cleanup system in 102 lb/hr sidestream of primary loop removes carbon and hydrogen 
impurities and removes or delays long-life volatile fission products in the helium by means of fil- 
ters, absorbers, gas chromatography, oxidizing beds, molecular sieves, and delay beds.'° 
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This decay heat exchanger protects both the 
downstream molecular sieve and the delay beds 
from localized overheating. 

The molecular sieves reduce the concentra- 
tion of CO, and water in the helium to less than 
1 ppm. After leaving the sieve beds, the 100°F 
helium is refrigerated to —-20°F, the operating 
temperature of the carbon beds. These carbon 
delay beds, which are packed with activated 
carbon and are maintained at —-20°F by refrig- 
erated ethylene glycol in external jackets, delay 
krypton for 43.6 hr (10 half-lives of ®Kr) and 
xenon for 1790 hr (14 half-lives of '**Xe). When 
the helium leaves the delay beds, the blower 
preheater heats it to 100°F and then the cleanup- 
system blower returns it to the primary loop. 


Critical Experiment 


A critical experiment used an assembly of the 
actual core, reflector, and fuel elements and a 
close simulation of the reactor control rods. 
First criticality with the assembly was achieved 
on June 9, 1965. 

The observed critical mass of 6.5 kg was con- 
siderably larger than the calculated 4-kg criti- 
cal mass." Studies of these results attribute the 
deviation from calculations to (1) neglect of the 
effect of streaming through the exit coolant hole 
and the fuel-element slot, (2) uncertainties in 
the total effective cross section of the fabricated 
graphite, (3) a need for more careful treatment 
of the neutron-thermalization and the self- 
shielding effects of the fuel elements, and (4) the 
absence in the critical experiment of the addi- 
tional reflector provided by the carbon insula- 
tion in UHTREX. 


On-Line Computer 


An on-line process computer will aid in the 
operation of the reactor, the reactor-coolant 
loops, and the coolant-purification system. 

The computer system includes a central 
processor with 4000 words of magnetic-core 
memory, 32,000 words of drum memory, an 
input — output typewriter, and apaper-tape punch 
and reader; there are four logging typewriters 
and one visual-display unit. The system sam- 
ples 600 analog and 300 digital inputs and sets 
24 analog and 444 digital outputs. Twenty 
interrupt channels allow time-sharing of input— 
output equipment and initiate special programs 
that depend on UHTREX system conditions.® 
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The computer samples analog channels ata 
rate of 200 channels/sec and monitors them for 
out-of-limit conditions. It can prepare logs, 
perform on-line calculations, keep reactor fuel 
and other system histories, and monitor reactor- 
operator activities. The computer also has the 
capability of sequencing valves in operations 
such as fuel loading and can control coolant flow 
and thermal output of the cooling system. Inad- 
dition, system diagnostic programs and emer- 
gency-shutdown programs will be available to 
assist in operation of the system.® 
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